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CHAPTER  I 


INTRODUCTION 


1 . 1 Background 

The  analysis  of  most  marine  structures  has  been  known  to  be 
extremely  complicated  and  difficult.  Thus,  it  becomes  essential  for 
one  to  idealize  marine  structures  by  simple  configurations  so  that 
an  approximate  and  yet  reasonable  estimate  of  a ship's  vibration  and 
acoustic  radiation  can  be  made.  For  instance,  the  motion  of  a 
ship  hull  reinforced  by  the  framing  elements  can  be  qualitatively 
treated  as  that  of  a plate  or  a shell  reinforced  by  beam  or  ring 
elements.  Likewise,  a vibrating  machine  mounted  on  the  beam  of  a 
ship  hull  can  be  idealized  as  a time-harmonic  point  force  excitation 
on  a beam-reinforced  plate. 

Vibrating  machines  have  been  known  to  be  supported  by 
foundations  attached  to  the  reinforcing  beams  of  a ship  hull.  This 
is  obviously  dictated  by  the  necessity  of  providing  a good 
structurally  strong  member  to  distribute  the  static  load  of  the 
machine.  This  procedure  of  attaching  machinery  to  reinforcing 
beams  has  proven  to  be  an  effective  method  for  reducing  the  hull 
vibration  and  the  resulting  radiated  noise.  This  is  due  to  the 
fact  that  the  impedance  termination,  as  seen  by  the  machine,  is  higher 
when  it  is  attached  to  reinforcing  beam  than  when  attached  to  the 
softer  ship  plating.  The  reduction  in  the  radiated  noise  when  a 
vibrating  machine  is  attached  to  the  reinforcing  beam  rather  than  the 
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plate  has  not  been  analytically  proven,  though  it  has  been  intuitively 
known  and  shown  in  some  experimental  studies. 

In  order  to  achieve  a reduction  in  the  radiated  noise  from 
reinforced  plates  and  shells,  it  is  proposed  to  attach  suitably 
chosen  and  judiciously  located  discontinuities  in  the  proximity  of 
reinforcement,  either  to  effect  sources'  mutual  destructive  interference 
or  to  change  the  acoustic  character  of  the  sources.  This  new  approach 
will  be  referred  to  as  "Acoustic  Fairing."  It  is  anticipated  that 
the  radiated  noise  from  a structure  with  a structural  discontinuity 
would  be  influenced  by  the  application  of  acoustic  fairing. 

To  effect  a reduction  of  the  radiated  noise  from  structures 
with  structural  discontinuities  by  the  application  of  acoustic 
fairing > one  must  first  obtain  a better  understanding  of  the  role  a 
structural  discontinuity  plays  in  the  radiated  noise  from  reinforced 
structures.  Thus,  one  must  obtain  estimates  of  the  radiated  noise 
from  a representative  structure  with  a typical  structural  discontinuity 
such  as  a plate  reinforced  by  a beam.  And  then,  the  method  of 
"Acoustic  Fairing"  can  be  applied  to  that  structure. 

1 . 2 Statement  of  the  Problem 

For  the  exploration  of  the  essential  physical  features  of  the 
"Acoustic  Fairing"  studies,  two  specific,  typical  problems  are 
investigated  within  this  thesis. 

The  first  problem  is  to  investigate  the  radiated  noise  from  an 
infinite  elastic  plate  reinforced  by  an  infinite  beam  excited  by  a 
time-harmonic  point  force  located  on  the  beam.  An  attached  beam  can 
be  considered  to  act  as  a structural  discontinuity  of  the 
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following  type: 

(1)  a mass  discontinuity,  which  will  normally  act  as  a line 
mass  impedance, 

(2)  an  elastic  discontinuity  which  acts  as  a path  of  energy 
propagation,  and 

(3)  a beam-like  discontinuity  which  acts  as  an  elastic  as  well 
as  mass  discontinuity. 

The  theoretical  derivation  and  numerical  computations  of  this  problem 
provide  the  background  material  for  the  assessment  of  radiated  noise 
due  to  a coupled  beam-plate  system. 

The  second  problem  is  to  devise  a simple  mathematical  model 
of  acoustic  fairing  by  attaching  an  additional  infinitely  long,  line 
spring-mass  system  to  the  left  and  the  right  of  the  beam  and  to  study 
their  interaction  with  the  existing  attached  beam.  The  radiated 
sound  pressure  and  power  from  a beam-reinforced  plate  are  computed 
for  various  values  of  the  sprung  mass,  the  spring  stiffness  and  the 
distance  of  the  systems  with  respect  to  the  attached  beam.  The 
analytic  solution  for  this  problem  together  with  the  numerical 
computations  yields  an  understanding  of  the  role  the  physical 
parameters  play  in  a successful  application  of  acoustic  fairing. 

1.3  Basic  Approach 

The  derivation  of  a single  governing  differential  equation  for 
the  basic  beam- reinforced  plate  is  partially  based  on  the  works  of 
G.  Lamb,  Jr.  (11  and  E.  E.  Ungar  [2].  It  begins  with  writing, 
separately,  the  equations  representing  the  flexural  displacements  of  a 
beam  and  a thin  plate  from  their  equilibrium  positions.  The  fluid 
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pressure  acting  on  the  beam-reinforced  plate  is  expressed  in  terms 
of  a velocity  potential.  By  employing  Fourier  transform  techniques, 
taking  advantage  of  the  symmetry  property  and  by  making  the  reasonable 
assumption  of  a narrow  attached  beam,  a single  coupled  governing 
differential  equation  is  finally  obtained.  The  integral  representations 
of  the  sound  pressure  and  power  are  set  up  by  substituting  the 
continuity  equation  at  the  interface  of  fluid  and  plate  into  the 
governing  differential  equation.  Again,  integral  transform  techniques 
are  employed  on  the  governing  differential  equation.  The  far-field 
acoustic  pressure  is  obtained  with  the  help  of  asymptotic  techniques. 
Numerical  solutions  for  the  radiated  pressure  and  power  are  made  with 
the  aid  of  digital  computer  IBM  370. 

1.4  Review  of  Previous  Investigations 

The  studies  of  vibration  and  acoustic  radiation  properties 
of  isotropic,  homogeneous  structures  such  as  membranes,  plates,  and 
shells  have  been  extensively  carried  out  in  the  field  of  structural 
acoustics.  Many  publications,  related  to  this  area  of  structural 
acoustics,  can  be  easily  found  in  the  literature.  Thus,  the  review 
of  this  section  is  confined  only  to  those  papers  dealing  with  the 
effects  of  discontinuities  or  non-uniformities  on  the  vibration  of  and 
acoustic  radiation  from  structures. 

G.  L.  Lamb,  Jr.  [1]  has  analyzed  the  vibration  of  an  infinite 
elastic  beam  attached  to  an  acoustically  unloaded  infinite  elastic 
plate  which  is  excited  by  an  harmonic  point  force.  He  expressed  the 
input  impedance  of  the  beam  reinforced  plate  in  terms  of  degree  of 
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coupling  of  the  beam  to  the  plate  and  the  ratio  of  wavelength  in  the 
plate  to  the  wavelength  in  the  bar. 

E.  E.  Ungar  [2]  has  studied  the  transmission  of  plate  flexural 
waves  across  reinforcing  beams.  He  developed  the  equations  for  the 
beam-plate  interaction  in  vacuo  and  expressed  the  plate  motion  in 
terms  of  reflection,  transmission,  and  near-field  coefficients. 

In  R.  Lyon's  work  [3]  on  the  sound  radiation  from  a beam- 
reinforced  infinite,  thin  plate,  he  computed  the  acoustic  radiation 
resistance  caused  by  an  incident  flexural  wave  in  the  plate.  He 
used  Ungar 's  results  [2]  for  the  reflected  and  transmitted  amplitude 
for  a plate  in  vacuo.  Thus,  the  influence  of  the  near-field 
acoustic  pressure  on  the  reflectivity  and  transmissibility  of  the 
attached  beam  was  totally  neglected. 

E.  L.  Shenderov  [4]  has  obtained  the  transmission  coefficient 
of  sound  through  a thin  plate  with  interjacent  supports.  He  found 
that,  for  certain  angles  of  incidence  of  a sound  wave  on  the  plate, 
a substantial  decrease  of  the  transmission  coefficient  will  result 
and  that  the  acoustic  transmissibility  of  structures  is  smaller  at 
low  frequencies  than  at  high,  which  is  valid  because  the  supports  are 
rigid.  D.  D.  Plakhov  [5]  has  considered  a similar  problem;  namely, 
the  transmission  of  a plane  sound  wave  through  an  infinite  laminated 
ribbed  plate.  He  claimed  that  when  the  mass  of  the  ribs  is  increased 
to  infinity,  the  results  of  calculating  the  transmission  coefficient 
agree  with  the  results  published  by  E.  L.  Shenderov. 

1.  I.  Klyukin  and  Y.  D.  Sergeev  [6]  have  studied  the  problem  of 


a flexural  wave  in  a plate  without  fluid  loading  scattered  by  point 
antivibration  systems.  They  showed  that  the  radius  of  influence  of 
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the  antivibrator,  inside  which  the  attenuation  is  greater  than  3 dB, 
is  about  0.3  to  0.4  times  the  flexural  wavelength  in  the  plate  for 
antivibrators . 

V.  T.  Lyapunov  [7]  has  discussed  the  wave  propagation  in  a 
liquid-loaded  plate  with  an  obstruction.  He  has  shown  in  this  work 
that  the  energy  transfer  through  the  liquid  depends  on  the  angle 
of  incidence  of  the  wave  and  the  ratio  of  its  frequency  to  the 
coincidence  frequency. 

I.  P.  Konovalguk  [8]  has  investigated  the  reflection  of  a plane 
sound  wave  from  an  elastic  plate  reinforced  with  many  infinitely  long 
stiffness  members.  He  mentioned  that  the  dependence  of  the  potential 
on  the  impedances  of  each  rib,  in  transverse  and  rotational 
vibration  mode,  has  a distinct  resonance  character;  the  sharpness 
of  the  resonance  depends  on  the  thickness  of  the  plate.  Also,  the 
flexural  impedance  of  each  member  plays  an  important  part  in  the  total 
reflection  of  sound  from  the  plate,  especially  for  large  angles  of 
incidence  on  the  plate. 

V.  N.  Romanov  [9]  has  analyzed  the  influence  of  reinforcing 
beams  on  the  radiation  of  sound  by  an  infinite  plate  driven  by  a 
field  of  random  line  forces  between  the  beams.  He  showed  that  the 
pressure  of  the  beams  yields  a considerable  increase  in  the 
radiation  at  low  frequencies. 

G.  Maidanik  and  E.  M.  Kerwin,  Jr.  [10]  have  considered  the 
acoustic  radiation  from  an  infinite  orthotropic  plate  and  from  the 
multipanel-structure,  including  fluid-loading  effects  and  structural 
damping.  They  pointed  out  that  heavy  fluid  loading  reduces  the 
acoustic  radiation  from  a point-  or  line-driven  orthotropic  plate. 
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It  also  shifts  the  resonance  frequencies  of  the  panels  in  the 
multipanel  structures  downward.  Structural  damping,  if  it  exceeds 
radiation  damping,  can  reduce  acoustic  radiation  from  free-wave  and 
panel-resonance  responses.  In  their  considerations,  they  dealt  with 
average  values  only,  for  instance,  under  the  assumption  that  when  a 
typical  spacing  of  the  ribs  is  smaller  than  the  flexural  wavelength 
in  the  plate,  the  ribbed  plate  is  considered  as  a homogeneous, 
orthotropic  plate.  And  when  they  discussed  the  radiation  efficiency 
of  the  panels,  a statistical  approach  was  adopted.  Therefore,  the 
specific  effects  of  the  individual  ribs  on  the  sound  radiation  was 
not  fully  understood. 

E.  5.  Jarmul  [11]  has  studied  the  sound  radiation  from  simple 
and  ribbed  plates.  He  experimentally  verified  that  a structural 
discontinuity,  such  as  a stiffening  rib  in  a finite  plate,  increases 
the  sound  radiation  from  such  a plate  when  excited  by  a randomly  time 
varying  point  force.  He  also  indicated  that  one  component  of  the  sound 
radiation  of  a finite  plate  can  be  attributed  to  a Helmholtz  resonator 
effect.  This  component,  as  D.  S.  Pallett  [12]  put  it,  may  arise  under 
certain  circumstances  because  of  the  interaction  of  the  force 
component  and  the  edge  component. 

M.  S.  Howe  and  M.  Heckl  [13]  have  examined  the  sound  radiation 
from  infinite  plates  due  to  the  interaction  of  bending  waves  with 
density  and  stiffness  fluctuations  in  the  material  of  the  plate.  They 
regarded  the  variation  in  density  and  bending  stiffness  as  random 
functions  of  position  on  the  plate.  They  also  assumed  that  for  low 
frequency  excitation,  i.e.,  below  the  coincidence  frequency,  the 
random  density  and  bending  stiffness  constitute  a stationary  random 
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process.  They  concluded  that  the  radiation  losses  brought  about  by 
the  phenomenon  of  scattering  due  to  density  and  bending-stiffness 
fluctuations  are  relatively  very  small  at  low  frequencies,  but  are 
comparable  to  loss  factors  caused  by  other  mechanisms  such  as 
internal  damping  and  energy  transfer  into  supporting  structures  at 
high  frequencies. 


CHAPTER  II 


MATHEMATICAL  MODEL  OF  ACOUSTIC  FAIRING 
2 . 1 Assumptions 

The  basic,  specific  structure  considered  in  this  thesis  is 
an  infinite  fluid  loaded  plate  reinforced  with  an  infinite  beam 
excited  by  a point  force.  Of  course,  it  is  known  that  infinite 
plate-beam  systems  do  not  exist  in  reality.  However,  the  justification 
for  considering  the  plate  and  the  beam  to  be  of  infinite  extent 
lies  in  the  assumption  that  the  structural  damping  of  the  beam- 
reinforced  plate  Is  sufficient  to  render  boundary-reflected  waves 
negligible  compared  to  the  waves  propagating  from  the  point  of 
excitation. 

Throughout  this  study,  it  will  be  assumed  that  the  plate  and 
the  beam  are  thin,  homogeneous  and  isotropic,  and  obey  linear 
constitutive  equations  [14].  For  the  structure  (plate  and  beam), 
stress  is  a linear  function  of  the  strain,  which  will  be  considered 
much  less  than  unity,  and  for  the  fluid,  pressure  deviations  from 
the  hydrostatic  pressure  are  a linear  function  of  the  density 
change,  and  thus  acoustic  pressure  is  just  a small  variation  from 
the  hydrostatic  pressure. 

Since  acoustic  motions  in  the  inviscid  fluid  are  always 
irrotational , e^en  if  viscosity  effects  are  considered  [15],  the 
problem  of  wave  mo" ion  in  the  fluid  thus  can  be  solved  by  a velocity 
potential . 

At  the  interface  of  the  structure  and  fluid,  the  fundamental 
boundary  condition  is  given  by  equating  the  normal  velocity  of  the 
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plate  to  that  of  the  fluid.  No  conditions  are  enforced  on  the 
tangential  velocity  component  in  front  of  the  thin  viscous  layer  at 
t'nt  boundaries  for  the  case  of  small-amplitude  periodic  motion,  this 
thin  ayer  acts  like  an  ideal  lubricant  [16]. 

In  most  physical  problems  posed  properly  by  a set  of  given 
boundary  conditions,  there  would  exist  only  one  physical  solution. 
However,  in  some  cases,  so-called  "paradoxes"  occur;  i.e.,  it  may 
appear  that  there  exists  more  than  one  solution.  For  instance,  the 
path  of  integration  in  the  integral  representation  of  a solution 
may  pass  through  a pole  of  the  integrand  on  the  real  axis,  resulting 
in  an  unbounded  value  for  the  integral  In  the  strict  Riemanian 
sense.  Although  the  integral  may  exist  in  Cauchy-Principle-value 
sense,  it  does  not  necessarily  represent  the  solution  of  a physical 
problem.  This  situation  is  often  a consequence  of  unrealistic 
assumptions,  usually  made  when  the  complete  equations  have  been  found 
too  difficult  to  solve,  and  an  approximate .oversimplified  set  of 
equations  has  been  substituted.  In  dealing  with  such  difficulties, 
the  physical  system  is  often  assumed  to  be  damped,  and  then  the 
solution  is  obtained  by  taking  a limiting  value  of  zero  damping. 

In  this  thesis,  the  classical  plate  equation  will  be  used  to 
develop  the  governing  differential  equations  of  the  point  force 
excited  beam-plate  system.  As  is  well  known,  the  classical  plate 
equation  does  not  account  for  shear  and  rotary  inertia  effects. 

Though  it  may  not  be  accurate  for  the  high  frequencies,  it  nevertheless 
will  adequately  illustrate  the  analytical  formulation  and  the  physical 
phenomena  [17].  The  width  of  the  attached  beam  and/or  distributed 
masses  will  be  assumed  smaller  than  the  flexural  wave  length  so  that 
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they  can  be  considered  as  mathematically  concentrated  stiffness  and 
masses.  The  formulated  governing  differential  equations  for  the 
vibration  of  and  acoustic  radiation  from  a beam-reinforced  plate 
excited  by  a point  force  are  those  involving  the  singularity  functions. 
The  desired  generalized  solutions  will  be  obtained  by  the  use  of 
integral  transform  techniques.  These  solutions  are  valid  for 
distributed  masses  and  stiffnesses , distributed  over  a vanishingly 
small  width,  but  can  be  considered  as  good  approximate  solutions  to  a 
thin  beam  attached  to  a plate. 

2 . 2 Classical  and  Generalized  Solutions 

The  solutions  of  boundary  value  problems  which  are  assumed  to 
be  sufficiently  smooth  and  to  satisfy  the  equation  at  each  point  inside 
the  region  of  definitions  of  these  problems  are  termed  as  "classical 
solutions"  [18].  Physical  phenomena  involving  discontinuities  which 
characterize  the  intensity  of  the  external  perturbation  can  be 
treated  conveniently  and  directly  by  the  use  of  the  singularity 
functions.  In  order  to  formulate  these  problems,  one  must  not  always 
insist  on  the  smoothness  of  the  solution  inside  the  region.  Thus, 
it  is  necessary  to  introduce  the  "generalized  functions"  to  generalize 
the  idea  of  a derivative  and  idea  of  a function.  The  solutions  of 
the  physical  problems  involving  singularity  functions  obtained 
by  the  use  of  generalized  functions  (or  distributions)  are  called 
'generalized  solutions." 

The  foundations  of  the  mathematical  theory  of  generalized 
functions  were  laid  by  Dirac  [19]  and  developed  by  Laurent  Schwartz 
[20],  Viewing  the  narrow  application  for  the  class  of  analytic 
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functions,  Schwartz  has  extended  the  concept  of  classical  functions 
and  proposed  a method  called  the  "Theory  of  Distributions"  which 
justifies  not  only  the  use  of  the  6-function,  but  also  the  use  of  all 
derivatives  of  the  6-function  [21].  He  considered  a generalized 
function  (or  a distribution)  as  a continuous  linear  functional  on  the 
space  of  testing  functions  which  are  continuous,  have  continuous 
derivatives  of  all  orders,  and  vanish  identically  outside  some  finite 
interval.  In  order  to  clarify  this  thought,  one  needs  to  define 
the  following: 

(i)  F(4>)  is  said  to  be  a linear  functional  on  the  space  of 
testing  functions  if,  to  every  testing  function  <]>(x),  a real  or 
complex  number  F(<£)  is  assigned  such  that 

F((f1  + <]>2)  = F (4^)  + F(<(>2) 

and 

F(A4>)  = AF($)  , 

where  A is  any  complex  number. 

(ii)  A functional  F(()))  is  said  to  be  continuous  if  the  sequence 
of  numbers  F (4:n)  converges  to  zero,  whenever  the  sequence  of  testing 
functions  <J>n(x)  converges  to  zero. 

The  following  are  some  important  properties  of  a specific 
family  of  singularity  functions  [22]  which  are  conventionally 
expressed  as: 
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Hn+1(x)  = 6n(x)  =0  if  x 4 0 


(2.2) 


and 


uu 

/' 


6“(x)  f(x)  dx  - (-l)n  f(n)(0)  , n-0,1.2,...,  (2.3) 


where  f(x)  is  a function  possessing  the  requisite  number  of 
derivatives . 

Some  useful  properties  of  6-function: 

(1)  A delta  (6)  function  acts  as  a derivative  of  a unit  step  function 


<5  (x)  = { H (x) } 


(2.  A) 


(2)  Sifting  property 


w 

I 


f(x)  6(x-x  ) dx  = f (x  ) (2.5) 

o o 


This  is  the  most  important  property  of  the  6-function.  Note  that 

this  property  is  satisfied  for  all  f(x)  independent  of  the  fact  that 

f(x)  has  to  be  a testing  function.  However,  f (x)  has  to  be 
differentiable . 

(3)  Equivalence  property 

Let  A^(x)  and  A^(x)  be  two  expressions  involving  6-functions 
and  ordinary  space  functions;  then  by  definition. 


A^x)  = Aj  (x) 


(2.6) 


if  and  only  if  J f(x)  A^(x)  dx  = J f(x)  A2(x)  dx  for  all  space 

— oo  — oo 

functions  f(x)  for  which  the  integrals  exist. 

(A)  Even  and  odd  property 


6n(x)  - (-1)"  6n(-x)  . 


(2.7) 


r 


A. 
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(5)  Higher  order  6-function  sifting  property 

OO 

J 6n(x)  f(x)  dx  = (-1)°  fn(x)|x=Q  = (-1)°  fn(0) 

—OO 


(2.8) 


To  summarize  this  section,  generalized  functions  can  be  handled 
algebraically  as  if  they  were  ordinary  functions,  provided  that  their 
use  can  be  justified  by  recourse  to  their  functional  representation. 
Also,  it  is  to  be  understood  that  if  an  equation  involving  generalized 
functions  is  multiplied  by  an  arbitrary  continuous  function  f(x), 
which  has  a sufficient  number  of  derivatives,  and  then  integrated 
from  to  with  the  functional  property  of  the  generalized  functions 
used  to  evaluate  the  integrals,  the  result  will  be  a correct  equation 
involving  ordinary  functions. 


2 . 3 Derivation  of  Governing  Differential  Equation  tor  a Beam-Plate 

System 

Consider  an  infinite  thin  elastic  plate.  Figure  (2.1), with  a 
thickness  h and  density  p lying  in  the  plane  z=0.  The  top  half-space 
of  the  plate  z < 0 is  a vacuum,  while  the  bottom  half-space  of  the 
plate  z > 0 is  occupied  by  an  acoustic  fluid  of  density  , sound 
velocity  c.  An  infinitely  long  beam  is  attached  on  the  top  of  the 
plate  at  x=0.  A time-harmonic  point  force  F^e  luJt  6(y)  is  excited 
on  the  top  of  the  beam,  where  F^  is  the  amplitude  of  the  force,  6(y) 
is  the  Dirac  delta  function,  and  (*)  is  the  circular  frequency  of  the 
harmonic  excitation. 

Starting  with  the  wave  equation  for  the  fluid  velocity  potential, 

♦ . 


(2.9) 
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letting  0 * 0e  " .and  using  the  rectangular  coordinate  (x,  y,  z). 
Equation  (2.9)  can  be  written  as  follows: 


| * k . o 

9x2  9y2  9z2 


(2.0) 


where  k^  = oj/c  (acoustic  wave  number). 


The  pressure  can  be  obtained  from  the  fluid  velocity  potential 


p(x,  y,  z,  t) 


3 

o 9t 


- icop  e 
o 


-iu>t 


(2.11) 


Denoting  the  transverse  displacements  of  the  reinforcing  beam 
and  the  plate  by  U(y,  t)  and  W(x,  y,  t),  the  coupled  differential 
equation  of  motion  for  the  beam  becomes  [1]: 


„T  9AU  92U  n 

EIb  TT  + \ 7T  " - D 
3y  9t 


— V2  W 
9x 


x=0 


- — V2  W 
+ 9x 


x=0 


+ Fq  e 1U)t6(y)  - bp. 


(2.12) 


3 

where  1^  = moment  of  inertia  of  the  beam  = 1/12  bh^  , see  Figure  2.1, 
* total  mass  of  the  beam  * pbh^  , 

D,  E,  Fq  and  p were  given  earlier.  The  bracketed  terms  in  Equation 
(2.12)  represent  the  shear  force  applied  by  the  plate  on  the  beam. 

The  equation  of  motion  of  a plate  coupled  to  an  acoustic 
medium  can  be  written  as  follows: 


2 -ft*  + 

3x2  3y2 


P 


< x < 00 


for , 


» < x < 


b 

2 


(2.13a) 

(2.13b) 


Let  us  define  the  Fourier  cosine  transform  pairs  in  one  dimension  x by 
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f ( 


(k')  = f 

‘ * Jo 

OO 

x)  '*/ 

^ A 


f(x)  cos  k'x  dx 
x 


f (k')  cos  k'x  dk* 

C X XX 


(2.14) 


-ik  x 
x 


and  the  Fourier  transform  pairs  for  one  and  two  dimensions  by 

O 

x JE 

f(x)  - — / f(k  ) 

SE  J x 


?(k  ) = — f f(x)  e A dx 

?E  J 


ik  x 
x .. 

e dk  , 
x 


(2.15) 


and 


f (k  , k ) 
* y 


2tt  / 


-ik  x -ik  y 

X V 

y)  e e ' dx  dy 


—•OO  — OO 
OO  OO 


1 f r~  ik  x ik  y 

f(x’  y)  " E J J f<kx’  ky)  e x e y dkv  dk,  . (2.16) 


x y 


— OO  _00 


Differentiating  Equations  (2.12),  (2.13a)  and  (2.13b)  with 
respect  to  time  t and  removing  the  time  factor  e_ilxJt,  one  obtains  the 
flexural  velocities  of  the  beam  and  the  plate  as  follows: 


Eib  ft  - “V  ■ - d 


9y 


fa  „2  , 

9 r,2 

*1 

< 

< 

^ V| 

L 1 

x=0 

o 

H 

X 

- ico  F 6 ( y ) 
o 


+ b OJ  <J>(z=0) 


(2.17) 


and 


bi^+2_A  + A 

l9x  9x  9y  9y 


2 2 
us  mv  = lo  (}>(z=0)  , 


for 


< x < 00 


® < X < 


(2.18a) 

(2.18b) 
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The  continuity  equations  between  beam  and  plate  along  the 
line  of  their  juncture  are: 


and 


u(y)  = v(x=0 , y) 

(2.19) 

y)  0 

9x  „ " ° • 

(2.20) 

Equation  (2.20)  is  a consequence  of  the  symmetrical  property 
of  the  beam-plate  configuration  and  the  drive  force. 

Applying  Fourier  cosine  transform  for  Equation  (2.18a)  on  x 

, u b 

for  the  range  'j  < x < 00 , one  obtains: 


k'4  v (k\  y)  + k’2 


X c X 


3 v(x,  y) 


x=0 


/ , d v (k' , y)  ~3  , . 

+ 2 -k*2  — l.v-(x».y> 


dy 


9x 


9x 


x=0 


x=0 

d*  vc(kx*  y) 


dy 


- U)2  m vc(k^,  y)  - Pq  O)2  <))c(k^>  y,  z = 0)  . 


(2.21) 


Rearranging  Equation  (2.21)  and  noting  that  d 


9x  9y 

which  can  be  deduced  from  the  Equation  (2.20),  one  obtains 


- 0 , 


x=0 


|>3y(x,  y) 

, 33v(x,  y) 

•I 

_ rv 

r 33v(x.  y) 

2 

9x  9y 

x=0 

, - D 
x=0+J 

a 3 

L dx 

x=0+- 

4 

* D ! k,H  v (k' , y)  - 2 k 

IX  ex7  X 


,2  ^ vc(k^»  y)  d*  vc(kx  y) 

o * 


dy 


dy 


2 2 
w m v (k',  y)  - p u)  <t>  (k ’ , y,  z * 0) 
c x 7 Ko  c x 7> 


(2.22) 


19 


v(x,  y)  cos  k’x  dx 
x 


D f—  V2  v(x,  y)|  ■ D k'^  f v(x,  y)  cos  k'x  dx 

9x  I x»o+  L x V 

OO  00 

,2  f 9"v(x,  y)  , , , . f 9^v(x,  y)  , , 

- 2k  I ir^~  cos  k x dx  + I J— ■ 1 — cos  k x dx 

x / -.2  x / . » 4 x 

V 9y  V 9y  J 

OO 

2 r 

- to  m / v(x,  y)  cos  k'x  dx 
0+ 

00 

2 r 

- Pq  oo  / <p(x,  y,  z = 0)  cos  k^x  dx  (2 


2 

- to  m 


y)  cos  k'x  dx 
x 


Similarly , 


- D f V2  v 
dx 


V2  v(x,  y) I = Dlk'^  / v(x,  y)  cos  k'  x dx 

x=o~  L x J x 

—.00 

® 2 ^4 

-2  k'2  / t-v,^x 2_.y cos  k'x  dx  + /"  9 cos  k'x  dx 

X <L  9y2  x J 9y4  x J 


- oo  m / v(x,  y)  cos  k^x  dx 


- p to 
o 


<p(x,  y,  z = 0)  cos  k^x  dx 


Substituting  Equations  (2.23)  and  (2.24)  into  Equation  (2.17),  one 
obtains  a universal  beam-plate  coupled  differential  equation: 


Dh1(/*/Jv(x- y) 


cos  k'x  dx\  - 2k'“ 
x / x 


-00  0 


[0  °°  J 

>2JI  r+  / |9  v<x>  y) 

x \/J  >2 

Noo  o+'  9? 


cos  k'x  dx 
x 


0 °o 

HI 


9 y(x,  y) 


cos  k 


x*  dx]  "U>2m  in  v(x,  y) 


cos  k'x  dx 
x 


i 


I 


' > 


IJ9 
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' I 

i 


o 00 


p OJ 

o 


| J + j |<J>(x,  y,  z = 0)  cos  k/ x dx  = - EIfa  - v^x  * ^ 
-oo  o+‘  * *y 


2 2 

+ co  ^ v ( x * 0,  y)  - iw  Fq  6 (y)  + b pQ  u <J>(x,  y,  z = 0) 


(2.25) 

Because  the  structure,  fluid  pressure,  and  excited  force  are 
symmetrical  about  x=0.  Equation  (2.25)  can  be  expressed  in  terms  of 
the  Fourier  transform  by  multiplying  both  sides  by  - — - . Moreover, 

SH 

due  to  the  fact  that  the  dimension  b is  a comparatively  small 
quantity  and  that  there  is  no  torsion  existing  along  the  beam,  one 
can  simplify  Equation  (2.25)  as  follows 


k ^ v(k  , y)  - 2k  2 

X x’  } X 


d2  v(k  , y)  d^  v(k  , y) 


dy 


dy 


- a)  m v(k  , y) 
x 1 


5o  ^ *(kx»  y»  2 = o)  = — f-  ei^-J^l — + w2Mv(x  = o,  y) 

JU  l 9y 


- iw  Fq  6 (y ) 


(2.26) 


where  El  * EI^  - bD,  M = - bm,  and  k^  is  replaced  by  k . 

It  can  be  easily  verified  that  Equation  (2.26)  is  the  Fourier 
transform  of  the  following  equation: 


J*v(x,  y)  . n9^v(x,  y) 

l 3«‘  8*23y2 


+ jV*.,  y> 

9y*‘  3y^ 


+ «(x)  ei  -a_vJjb  y) 
3y4 


- tv  [m  + M 6(x)]  v (x , y)  = - iw  F 6(x)  6(y) 

o 


+ Po  w 0(x,  y,  z = 0)  . 


(2.27) 


* ^ 
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Equation  (2.27)  represents  the  coupled  beam-plate  differential  equation 
which  will  be  used  to  develop  the  governing  differential  equation 
for  the  acoustic  fairing  problem  in  subsequent  sections. 

2 . 4 Far-Field  Pressure  of  the  Beam-Plate  System  with  Two 
Additional  Spring-Mass  Systems 

2.4.1  Statement  of  the  Problem.  By  providing  two  additional 
spring-mass  systems  symmetrically  about  x=0  as  shown  in  Figure  2.2,  the 


governing  differential  equation  for  the  whole  system  will  be,  of 
course,  different  from  that  derived  in  Section  2.3.  However,  the 
desired  differential  equation  for  this  simple  idealized  model  can  be 
obtained  by  modifying  the  previous  equation  without  resorting  to  the 
rederivation  of  the  new  coupled  plate-beam-sprung  masses  system. 

This  can  be  done  because  the  chosen  simple  model  for  the  sprung 
masses  does  not  involve  complicated  couplings  with  the  plate. 

The  following  differential  equations  represent  the  coupled 
motion  for  this  problem: 


(beam-plate):  ~4- 

D 6 W 


4 + 20  2 

9x  9x  9y 


2^5  + D + S(I)  EI 

3y4  Sy‘ 


+ [m  + M 6(x)]  — ^ = F e 6(x)  6(y) 

9t2  ° 

+ k^(W^-W)  6(x+d)  + k^Wj-W)  6(x-d)  - p(x,  y,  z = 0) . 

(2.28) 


(spring-mass): 


mx  W1+  k1(W1-W) 


0 , 


(2.29) 


where  represents  the  motion  of  the  mass  M^,  see  Figure  2.2,  k^ 
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the  spring  constant  and  d is  half  the  separation  distance  between  the 


masses . 


Combining  Equations  (2.28)  and  (2.29),  and  differentiating 
the  result  with  respect  to  time  t,  one  obtains  the  following 
differential  equation: 


J + 2 — | 55-  + + <5(x)  El  — ? - cj2[m  + M 6(x)]  v 

*\  ^ 


4 2 2 4 

9x  9y  9y 


2 k.  w 

= -iuj  6(x)  6(y)  + pQ  w (x , y,  z = 0)  + — ^ — 2~  v 6(x+d) 


U1  -u>  ■ 


k uj 

+ — 2 — 2 v ^x"d) 

00^  -W 


2.4.2  General  Solutions.  Applying  the  Fourier  transform  on 
variables  x and  y of  Equation  (2.30),  one  can  get 

D[(-ik  )4  v + 2(-ik  )2  (-ik  )2  v + (-ik  )4  v]  + — (-ik  )4  v(x=0 , k ) 


2 - oj2  M . 

uj  m v v(x=*0,  k ) 

/2?  y 


"2tt  + po  “ ♦ 


1 / k u>  1 ik  d 

+ — 5 — ~ v(x  = -d,  k ) e x 

^L12-U2I  y 


. k co  \ -ik  d 

+ - — ■ ■■  v ( x *5  d k ) e X 

2 2/  V^X  a’  V e 

/2tt  'u)^  -to  1 7 


(2.31) 


In  order  to  solve  the  complete  problem,  the  velocity  potential 
will  prove  to  be  a convenient  tool.  Starting  with  the  wave  equation 
on  the  velocity  potential  for  the  acoustic  medium,  given  in  Equation 
(2.10), 
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i±  + i±  + ii  + k 2 * , o 

3x2  3y2  9z2  ° 

and  applying  the  double  Fourier  transform  on  the  independent 
variables  x and  y of  Equation  (2.10),  one  obtains: 


(2.10) 


-4  + (k  2 - k 2 - k 2)<t>  = 0 . 
dz2  0 x y 


(2.32) 


The  solution  of  the  differential  Equation  (2.32)  for  outgoing  waves 
can  be  written  as: 

izX~* 


<t> (k  , k , z)  = A(k  , k ) e 
x y x y 


2 2 
k - k Z 
x y 


(2.33) 


For  the  determination  of  the  coefficient  A,  one  must  solve  for 
d4>(k  , k , z = 0) 


X ^ , 4>(k  , k , z = 0),  v(x=0,  k ),  and  v(x  = +d,  k ) , 


dz  ’ T ' x ' y ' ' ' y 

and  then  substitute  these  into  Equation  (2.31).  Thus,  the 

transformed  potential  at  the  surface  of  the  plate  is  given  by: 


<{>(k  , k , z = 0)  = A(k  , k ) 
x y x y 


(2.34) 


Using  the  continuity  condition  of  the  normal  velocity  at  the 
plate-fluid  interface. 


v(x>  y,  * - 0)  . -M* 


(2.35) 


z=0 


the  following  transformed  equations  results: 


v(k  , k , z * 0) 

x y 


d<J>(k  , k z) 

* y - 


dz 


I z*0 


/ 2,  27  2 

r— o o o 

-Ai/t  2-k  2-k  2 e ° X y 
o x y 


= -Aiv'k  2-k  2-k  2, 
z*0  o x y 


(2.36) 
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v(x,  k , z = 0)  = 

/2n 


OO  — 

— f Ai/k 

2?  J 

J—co 


2 2 2 lksx 

-k  -k  e dk 
s y s’ 


(2.37) 


v(x  = 0 , k , z * 0) 

y /2tT 


= - — f A(s,  k )v/k  2-s 

m J y ° 


2,2, 

s -k  ds 

y 


(2.38) 


and 


v(x  = + d,  k , z = 0) 

/2tt 


= — f A(s , k )i/k  2-k  2- 

y o y 


s2  e±lsd  ds  . (2.39) 


Substituting  Equations  (2.35)  and  (2.39)  into  Equation  (2,31),  one 
obtains : 

4 °° 


^(k  4+2k  2k  2-t-k  4)  - 

ILillt 

y 

f A(s , 

2tt 

OO 

x x y y 

n 

i 

A 2-k  2-k  2 ♦ “ ,M1 

O X y ^TT 

f A(s  , 

k )Jk  2 
y o 

2 2.2, 

-s  -k  ds 

J y 


-iw  F 


a o 2.  l|  1 ^ x f . , , , -isd  / 2 FT 

■2F-  + p0  wA-d— r~r  6 A(s’V  6 Vko  _s  "N 

\ to,  — CO  » ' J 


ds 


,|  k,  U)  \ -ik  d /*  , < — 

Trr|  2 2 ) 6 X /A(s>  V el"V^ 

>oj.  -w  ' » 7 


2 2 2 

s -k  ds  . (2.40a) 

y 


Rearranging  Equation  (2.40),  an  integral  equation  on  the  coefficient 
"A"  results: 


I iw  F 2 

. \ o , w Mi 

A(kx>  V “ [~2F“  + 2tt 

4 


ou 

f A(s,  k ) A 2-s2-k  2 ds 

J y o y 


i El  k 


s2-k  2 ds 

y 


i | 

k ur 

2tt| 

2 2 
-U) 

X—  f A(s , k ) i/k  2- 

J y o 

i4)«‘v  /*<••  y «'i8Vko2-*2- 


k ds 

y 
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i 

k u) 

2tt 

2 2 
'co,  -ui 

k ) elsd/k  2-s2-k  2 ds 

y V o y 


,n  / 2 . 2 . 2 L 2 ^ , 2,2  u)  ml  2 

lD/k  -k  -k  (k  + k ) — — + p u 

o x y x y D o 


(2.40b) 


Before  solving  this  integral  equation,  one  can  adapt  the 
following  notations  to  expedite  the  manipulation;  namely. 


Loo  F 2 ik,  ur  . 

o _ _ iw  M 1 i El 

2v  * 81  2tt  ’ g2  " , , 2 2 ’ Sy  ° 2ti  ’ 

2tt(co1  -w  ) J 


)v^  2-k  2-k  2 f(k  2 + k 2) 


Q(k  , k ) = iD/k  -k  -k  (k  + k ) 
x y oxy  x y 


+ p U) 
o 


(2.41) 


Now,  Equation  (2.40b)  can  be  rewritten  in  the  following  brief  form: 

, ik  d -ik  d 

' + *lCl  - Sy\  G1  + *2*  G2  + *2«  * G3 


A(k  , k ) = 
x y 


Q(k  , k ) 
x y 


,(2.42) 


(k  ) = f A(s , k ) 
1 y J y 

— oo 
00 

G2(V  = / A(s’  V 

—OO 
CO 

3 (ky)  = j A(s,  ky) 


. , -isd  f 2 2 7 2 , 

k ) e Ik  -s  -k  as 

y V ° y 


G_  (k  ) = / A(s , k ) elsd  AT2-s2-k  2 ds  . 
3 y / y Vo  y 


(2.43) 


It  is  expected  that  is  equal  to  C^,  but  for  the  time  being, 

G^  and  G^  are  to  be  considered  as  two  unrelated  unknowns. 

Substituting  the  expression  for  A in  Equation  (2.42)  into  the  integrals 
in  Equation  (2.43),  three  independent  equations  can  be  set  up  to 


solve  for  G , G0 , and  G.  as  follows: 


<f  * eici  - syky  ci>  J[  *♦  — ocJTTTr-2 

pe-is<i  fo2-s2-k  2 

+ g2C3  J ■ Q(s,°k  ) ^ = (f  + «1G1  - VyV1!  + 82G2I 

— ao 


+ g2C312 


(f  + g G - s k G 
11  y y 


° -isd/  22,2 
-e  _/k  -s  -k 


) fl— V! o. 

1 J Q(s , 


~ -isd  isd/  2 2,2 

re  e /k  -s  -k 

+ 82G2  / Q(s  , IT)  *“  ds 

— oo  y 


00  -isd  -isd  I.  2 2 , 2 

/e  e >k  -s  -k 

Q(57rr ^dS  ' “ + SlGl  - Syk/Gl)Ii 

— ao  y 


f 82G211  + g2G3 1 3 


(f  + glGl  - Syky4Gl)  J~ 


eisd/k  2-s2-k  2 
V o y 


Q(s,  ky) 


feisd  elsyko2-s2-k  2 } eisd 

' 82G2  J Q(s,  O + g2G3 / qTs,  k )' 

-°°  y y 


‘ (£  + glGl  - Syky  G1)I2  + h°2l3  + g2G3Il  * 


/k  -s  -k 
o \ 


Q(s,  k ) 

y 


isd  l 2 2 , 2 


00  -isd  A 2 2 . 2 


Q(s,  k ) 


Q(s,  k ) 
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(2.47) 

7~ 2 2 2 

Since  /k  -s  -k  and  Q(s,  k ) are  even  functions  with  respect 

° y y 


to  s,  one  can  prove  the  following  relationship: 


I3  = I3  . (2.48) 

Substituting  the  above  relations  into  Equations  (2.44),  (2.45), 
(2.46),  one  forms  the  following  three  independent  s.multaneous 
equations : 

(1  - gllx  + syky4I1)G1  - (g2I2)G2  - (g2I2)G3  = fl3  , (2.49a) 

(-g^2  + Syky4l2)Gl  + (1  " 82I1)G2  " (82I3)G3  = fI2  (2-^b) 

and 

(~g1I2  + syky  I2)G1  - (g2I3)G2  + (1  8213)63  = fl2  . (2.49c) 

By  inspecting  Equation  (2.49b)  and  Equation  (2.49c),  one  can 
easily  conclude  that 

G2  - G3  . (2.50) 

The  determination  of  the  coefficient  "A"  in  Equation  (2.42)  can  be 
accomplished  readily  if  one  takes  advantage  of  the  property  of 
equality  in  Equation  (2.50).  Thus,  when  Equation  (2.50)  is 
•ubstituted  back  to  Equation  (2.42),  the  following  simpler  relation 
results : 
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where 

ik  d -ik  d 

2g2G3  cos  kxd  = g2G2  e X + g.,G3  e X . 

Following  the  same  procedure  stated  above,  a new  set  of 
algebraic  equation  on  G3  and  G3  are  obtained: 


°/  2 2 , 
• /k  -s 


<:,-(£*  gjCj  - syky\>  J -°c-  ds 

— OO  y 

oo  / ' n ^ 2 

/cos  k d k “-s  -k 

— , r * , (f  + glCl  - . k/opij 


* 2‘2°3l2 


G3  • <£  + *l°l  ' J 


lsd  / 2 2,2 

e ^/k  -s  -k 

lx  o X_ 

Q(s,  ky) 


+ 2*2°3  J 

— oc 

+ 282G3I4 


isd  , / 2 2,2 

e cos  sd  /k  -s  -k 

~Y  ° V 

Q(s,  ky) 


ds  * (f  + g G, 


Syky  Gl)l2 


00  dsd  i t 2 2 , 2 

/e  cos  sd  Ik  -s  -k 

V ° v 

Q (s , k ) 

— oo  y 


/V*  - 

”qu7 


-s  -k 


° / 2 2 2 

cos  sd  /k  -s  -k 
Vo  y 


Q(s,  k ) 

y 


Vs:!/k  2-s2-k 2 
Vo  y 


Q(s,  k ) 

y 
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If  one  sets 


’/ 


00  2isd  / 2 2 , 2 

e /k  -s  -k 

JL° ds  , 


Q(s , k ) 

y 


(2.53d) 


one  can  express  the  integral  I.  in  terms  of  I and  I-  as  follows: 


^ =2(I1+I3)  * 


(2 . 53e) 


Rearranging  Equations  (2.52a)  and  (2.52b),  one  gets 


(1 


- Hh  * W VG1  - 2*2:2G3  ‘ fIl 


and 


(-  glX2  + Syky  VG1  + U - 2g2I4)G3  = £I 


2 • 


(2.53a) 


(2.54b) 


The  two  unknown  functions  and  in  the  above  two  simultaneous 
equations  can  be  easily  solved  by  using  Cramer's  rule. 

The  determinant  of  the  coefficients  is 

A = (1  - g1I1  + s k V)(l  - 2g2I4)  + 2g2I2(-  glI2  + s k 4I2)  . 


(2.55) 


Solving  Equation  (2.54)  for  the  two  unknown  functions  G^  and  G^, 
the  following  expressions  are  obtained: 


51  = A(I1  " 2gl1l14  + 2g2£2  > (2'56a) 


f 1, 


and 

G3  ■ f(I2  - *lllh  + Syky4llI2  + * (2.56b) 

Finally,  the  coefficient  "A"  is  obtained  in  the  following: 


A(kx’  ky>  * ^ +VIr2g21lV2g2I2  ) ’ SVWlVVl  1 


+ 2g2  I2  cos  kxd )/A] » Q^(1-g21i_g2I3+282I2  cos  kxd) 


(2.5?) 


31 


The  radiated  acoustic  pressure  can  be  obtained  from  the  fluid 
velocity  potential  by  the  relationship  in  Equation  (2.11),  one  gets 
an  integral  expression  as  follows: 


p(x,  y,  z,  t)  = -ioo  p <J>  e 


-iwt 


where 


/ 2 2 2 

iz/k  -k  -k  i(k  x+k  y) 


• r a,°°  o x y x y'  _ ,, 

iwp  f . ^ re  e C dk  dk 

o -iu)t  / / ms  x y 

2 IT 


‘//! 


“-“IdX.  2-k  2-k  2 [<k  2+k  2)2-  - 

o x y x y 


m-|  2 

D~  "PoU 


oo  oo 

,1  2 2 -iwt  ' r& 

= (2?}  PoW  Foe 


X2  • 2 ’ 2 


If 


lz/k  -k  -k  i(k  x+k  y) 

0 * y e * y C dk  dk 

ms  x y 


- 2,2,  2 , . 2 , 2 2 U;2m,  2 

•°°  -°°iD/k  -k  -k  [ (k  +k  ) — ]+p 

o x y x y D o 

(2.58) 


Cms  = A(1  " 82X1  ~ *2*3  + 2g212  COS  kxd)  ’ 


2.4.3  Evaluation  of  Far-Field  Pressure  by  Asymptotic  Methods. 

In  order  to  carry  out  the  integration  in  Equation  (2.58),  the  following 
independent  variable  transformations  are  made  [23]: 

/~~2  2 2 

k = k sin  cx  cos  8,  k = k sin  a sin  6,  k = /k  -k  -k  = k cos  a 
xo  yo  zoxyo 


x = R sin  4>  cos  0,  y = R sin  (p  sin  0,  z = R cos  <p. 
The  Jacobian  of  such  a transformation  becomes 


(2.59) 


3k 

3k 

X 

X 

9a 

38 

3k 

3k 



__Z 

3a 

36 

k sin  a cos  a. 
o 


(2.60) 
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Therefore , 

dk  dk  = k ^ sin  a cos  a da  d3  . (2.61) 

x y o 

The  integration  with  respect  to  3 will  be  performed  between 
0 and  2tt,  while  the  integration  with  respect  to  a cannot  be  restricted 
to  real  values  of  this  angle.  The  integration  limits  for  the  argument 
a can  be  determined  as  follows: 


(i) 

when 

k = k =0:  k = k ,a  = 0 

x y z o 

(ii) 

when 

k -*•  + <»  and/or  k -*•+“:  k ■»  l®  . 

(2.62) 

x — y — z 


since  cos  a = k /k  = i 00  , one  has  oi  = •?•  - i°°  . 

z o 2 


As  a result,  the  integral  in  Equation  (2.58)  transforms  to  the 
following  integral: 


lk  R cos  4>  cos  a *jr  iknR  sin  <J>  sin  a cos  (8-0) 
P = P k F / e ° {j 


. (1  )2  „ A f 

2tt  o o o J 


«/>' 


o 


V(a,  8)  dB)  sin  a da  , 


(2  63) 


where 

cos  a C (a,  6) 

V(°-  8>  ~ Q(».“) • (2-64) 

Since  the  exact  solution  for  the  radiated  pressure  integral 

representation  in  Equation  (2.63)  cannot  be  obtained,  an  approximate 

solution  is  desirable.  For  the  solution  of  the  far-field  pressure 

when  k R >>  1 and  k z >>  1,  a useful  mathematical  tool  called  the 
o o 

method  of  steepest  descent  will  be  used  here  [24,  25,  26].  The 
fundamental  concept  of  this  asymptotic  method  is  to  distort  the 
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original  integration  path  into  the  complex  path  of  steepest  descent 
(PSD)  so  that  the  complex  line  integral  can  be  reduced  to  evaluating 
a single-valued  function  at  the  saddle  point  where  it  receives  the 
major  contribution  in  the  whole  distorted  path  [27].  It  should  be 
noted  that  any  singularities  crossed  in  the  deformation  of  the 
original  integration  path  into  the  path  PSD  should  be  included  as 
discrete  pole  contributions  according  to  the  residue  theorem  [28]. 
However,  these  discrete  residue  terms,  in  general,  represent  the 
near-field  pressure. 

To  obtain  an  asymptotic  solution  by  this  method,  one  can  let 
g(B)  represent  the  part  of  the  power  of  the  exponential  term  in 
Equation  (2.63);  i.e.. 


g(B)  = i[sin  sin  a cos  (8-0)]  . (2.65) 

The  saddle  point  8q  is  the  root  of  the  following  equation: 

g ' (B  ) = - i sin  sin  ot  sin  (8  -0)  = 0 , or  B =6,  tt+0 
o o o 

(2.66) 

If  one  puts  B - B1  + i&2  » then, 

g(B)  * sin  <J>  sin  a[sin  (8^-0)  sinh  8^  + l cos  (B^-0)  cosh  3 • 

The  regions  of  convergence  of  the  integral  in  Equation  (2.63) 
are  determined  by  real  part  of  function  g(B) • namely, 

Re{g(B)}  ■ sin  <$>  sin  a sin  (Bj-0)  sinh  • (2.67) 

In  the  upper  half-plane  > 0) , the  regions  of  convergence 


are : 
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sin  (8^-0)  < 0 
or 

-7T+0<ei<0  , TT  + 0<6i<2Tr  + 0 

And  in  the  lower  half-plane  < 0) » the  regions  of  convergence 

are 

sin  (8^-0)  > 0 
or 

0 < [3  < it  + 0 

The  steepest  descent  paths,  which  are  found  from 

lm{g(3)}  = constant  (2.68) 

pass  through  the  two  saddle  points  within  the  region  between 
3 = - ~2  and  3 = are  then  defined  by  the  following  equations 
(see  Figure  2.3): 

(i)  For  saddle  point  at  (0,  0)  : cos(8^  - 6)  cosh  * 1 • 

(ii)  For  saddle  point  at  (0+r,  0)  : cos(B^-0)  cosh  82  * ~ 1 

The  asymptotes  of  steepest  descent  paths  are  determined  from 

cos(3j-0)  -*-0  , or31  = --|^-+0  , - y + 0 , j + 0 >‘f~+9  • 

Thus,  it  can  be  seen  that  the  suitable  steepest  descent  paths 
for  the  present  case  can  be  taken  from  (-  tJ-  + 0 + i00)  to  (0) , then  to 
("2  + 0 - i““) , and  then  from  (y  + 0 - i°°)  to  (tt+0)  and  then  to 
(■|—  + 0 + i00)  . In  other  words,  these  steepest  descent  paths  can 
be  shown  to  be  equivalent  to  the  inner  integral  in  Equation  (2.63), 
provided  that  there  are  no  singularities  within  the  complete  close 
contour  described  by  (-  ^ + 0)  > (0)  •+  (tt+0)  + 0)  ■*  (y-  + 0 + l00) 
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(-c2)  - (-c  ) - (-  j + 0 + i°°)  - (-  y + 0) . 

The  proof  of  this  statement  is  simple,  since  the  periodicity  of 


the  integrand  is  2 tt  ; thus  , 

t 

2tt  2n  0 

(1)  [l  ] dB  = [ + f 


j[  } d6  = If  + f - J l [ ] dB  = J 

0 0 *'ir  ^ 

_ _ + 6 . _ + 0 - j 


[ 1 d6 


-!+  e 


[ 1 dB- 


I ] dB  , 


- 2"  + 0 + i® 


+ 9 + i°° 


- 7 + 0 ~ + 9 + i00 

j[  ♦/ 

- y + 0 + i°°  ~+0 


M dB  = 0 


From  the  residue  theory. 


ff-e  ?+e+i» 


6 [ ] dB 


- ~ + 9 + 9 


-/  -/*  1 

C-  C.  -TT 


>[  ] dB  = 0 . 


C2  C1  -j  + 0 + i=°l 


Therefore , 


[ ] dB 


[ 1 dB 


■hi 


t 1 dB  , 


C1  c2 


where  [ ] denotes  the  integrand  in  Equation  (2.63). 


Letting 


u = g(3Q)  - g (B)  , 


u^  * i sin  <ji  sin  a[l  - cos(B-9)]  , on  path  (2.69a) 


u„  = i sin  <)>  sin  a[-  1 - cos(B-9)]  , on  path  C2,  (2.69b) 
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then,  the  inner  integral  in  Equation  (2.63)  can  be  written  as  follows: 

ik  R{i  sin  <p  sin  a - u 


TT  IT 

-j  + 6 - i°°  - — + 0 + i°°. 

7 -/ 

L 0 6 


on  C, 


ry-  + 0 + i°°  - y + 0 + i“ 


I -1 

TT+0  TT+0 


1 V(cx,  6)  ~ du. 

du^  1 


ik  R{-i  sin  <t  sin  a - u„)  ,Q 

2 »c*.  « 


on  C, 


(2. 70) 


where  and  are  two  paths  of  integration  as  shown  in  Figure  2.3. 

Before  the  integral  in  Equation  (2.70)  can  be  simplified  any 

further,  the  terms  4^—  and  have  to  be  expressed  in  terms  of  u,  and 
du^  du2  r 1 

u2 , respectively.  If  Equation  (2.69a)  is  differentiated  with 
respect  to  u^ , one  obtains 


du  i sin  $ sin  a sin(0-(}>) 

1 /2i  sin  <p  sin  a u^  - 

Since  k^R  is  assumed  to  be  large,  only  small  values  of  will 

contribute  significantly  to  the  integral.  Thus,  the  above  expression 
<J  ' 

for  yy—  can  be  expanded  into  a Taylor's  series  in  and  by  retaining 
the  leading  term  in  the  series,  an  approximate  expression  for  the 
derivative  results: 


dB  _ / 1 -1/2  -Tri/4  / 1 -1/2  3m/4 

du^  / 2 sin  sin  a U1  e ’ °r-y  2 sin  <p  sin  a U1 


For  u:\tqueness,  one  requires  that  the  above  expression  with 

— TTi/4  TT 

the  factor  e is  used  for  the  path  from  0=0  to  0 = y + 0 - i00  , 

while  that  with  the  factor  e^  used  for  the  path  from  0«0  to 

0 = - y + 0 + i°°. 
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Similarly , 


dB 

du„ 


-1/2  Tri/4 
e 


/ 1 

/2  sin  <j>  sin  a U2 

should  be  used  for  the  integration  path  from  rr+0  to  -jp-  + 0 + i®  and 

dB  s / 1 u- 

du^  ^1  2 sin  (p  sin  a 


-1/2  5tti/4 
e 


used  for  the  path  from  8 = tt+0  to  8 = + 0 - i°°  . 


Thus , 

2tt 

ik  R sin  <p  sin  a cos  (8-0) 


I 


V (a,  8)  d8 


ik  R sin  <t>  sin  a, 
2V(a,  0)  e ° 


7 1 1 -1JM  f -1/2  „ 

_7'  sin  1 sin  a / e U1  U1 

V 0 

1 7 r"  ~~  tti/4  f "k0Ru2 

_/2sin$sina  e J e u 


+ 2V(a,  tt+0)  e 


-ik  R sin  <J>  sin 


1/2  a 
2 du2 


Hi 

VkoR  SJ 


sin  <J> 


V(a,  0)  e ° 


ik  R sin  <t>  sin  a 


1 -TIi/4 

. e 

sin  a 


— -ik  R sin  4>  sin  a / 1 TTi /4 

V(a,  tt+0)  e ° /sin  a e 


(2.71) 


Substituting  Equation  (2.71)  into  Equation  (2.63),  an  integral 


representation  for  the  far-field  fluid  pressure  p is  obtained: 

n. 

— — X°° 

2 /'  ■ \ ■ 9 2f  ik  R cos  0 cos  a 

p . p JvTT^  K 1 e°  lv<“’e> 

2tt  "o  oy  o r J 


(a. 


ik  Rsin4>sina  » , ... 

•”i/  ♦ «•>  • ° 


ik  Rsin<f>sina 


1 Tti/4 
e 


sina 


sin  a da 


(2.72) 


N>|  ^ 
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-100 


/ 


Noting  that  V(a,  0)  = V(a,  tt+0)  and 


lk  Rcos4>cosa 


ik  Rsintfsina  / — — .. 

V( a,  tt+0)  e ° ~ysin  i ^ Sln  a 


o 

‘ l 


ik  Rcos<i>cosa  ik  Rsin<t>sina  r—. . .. 

e 0 V(a,  0)  e ° sin  a da  , 

’ ysina 


+ ioo 


the  expression  in  Equation  (2.^2)  can  be  rewritten  as 

— — loo 
2 

1.  »2  0 „ , 7- C ik  Rcos^cosa 

h)  vV'o-p?  I • ° «»■  6> 


- ~ + i» 


ik  Rsin<j)sina  _ . 

O -TTi/4  J— — , 

e e rsina  da 


1 _ ioo 

2 


N2„  2,  2 TTi/4  / f 

F0PoW  ko  e ykoRsin4)  J 


ik  Rcos(a-<))) 

> 0 V(a,0)v/s  ma  da. 


_ _ + i°o 

(2.73) 

The  integral  in  Equation  (2.73)  can  be  obtained  by  following 
tlie  same  procedure  and  method  outlined  above.  Let 

s(a)  = i cos(a-<}>)  , 

then,  the  saddle  point  ao  is  defined  by 


letting 


Then 


s'  (a  ) = - i sin(a  -(())  =0  , or  a = <J>  , 

o o o 


u~  = s(a  ) - s(a)  = i - i cos(a-<f>) 
3 o 


da  _ /r  -i/2  -ni/4  _ rr  -u 

du  V2  u3  e * ofV2  U3 


-1/2  3 rri/4 

e 


and 
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i(k  R-tt/4) 

, ° q\  "^i/A  2 

V(<J>,0)  e — /sin 

/2 


ik  F 
o o 


ik  R 

o , 

e cos  9 


2ttR  ik  mcostf) 

o r,  to  . 4. , 

~ n 11 9 sin  <f>l 


C (0,  0)  , 


(2 


1 - 


where 


C (<fr,0)  = 

ms 


ml  \l  “l 


0 „ 0 


1 + Pll  2 II  2~  ~2  K2I?~  cos  (k  dcosOsinif)  - I.  U-I  *" 
'p  h M co-  -w  » I -1 

o ] 


1 " Vi 


0)  M 


ml  1/  “l 


pS  looh2/l 

o o 1 


mi  1 

ui  | 

Ip  h2l 

2 2 | 
W.  -CO  1 

+ p,  U sin  4> 
3 y 


I,0'  * I °I  ° - 21  °21 
1 1 Z 2 


l_l ll 

1 "'2  1 

Ip  h2/ 

1 2 2 1 

1 00  —to  » 

“i-imi yV12^2)  ■ v-VT- 


u3  = 12HH(;H  u-v  > , 


o o 


Uy  ’ q3^1-v  )(q2  + 3<i2  + 3q2^ 

» . 

[ 0 _ J A.  - g 2 - n2  dn  


ia  v l-q 

2 

-I 

u 


b - ^ 


X 2~2 

/l-q  -n 
2 


{b(n2  + q2)  - 1}  + 1 


V2}/12 (1-V2)  , 


, q » sin  0 sin 


n‘ 
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■/ 


£_ 

Po 

ik  dn 
o 


sJ'V 


12p(l  - y2) 
E 


A 


2 2 . 

g - n dn 


ia/l 


q2  - n2  (b(n2  + q2)  - 1}  + 1 


and 


H 


2ik  dri  y- 
e ° A 


2 2 
q - n dn 


ia/1  - q2  - n"  ib(n2  + q2)  - 1}  + 1 

2.4.4  Discussion  of  the  Possibility  of  Replacing  the  Original 
Integration  Paths  by  Paths  of  Steepest  Descent.  To  deform 
the  original  integration  paths  into  paths  of  steepest  descent 
continuously,  certain  singularities  may  be  encountered  in  the  complex 
a and  complex  3 planes,  which  are  enumerated  below: 

(i)  /sin  oc  has  a branch  point  at  a = 0.  However,  this  point  is 
of  no  great  significance  since  the  steepest  descent  path  shall  not 
pass  through  it. 


(ii)  In  the  present  case,  the  only  singularities  present  would 
be  poles  of  the  function  V(a,  6)  in  the  6-plane  and  those  of  the 
function  V(u,  <p)  in  the  a-plane.  If  such  a pole  is  crossed  in  the 
deformation  of  original  integration  paths  into  steepest  descent 
paths,  then  the  residue  at  the  pole  must  be  added.  The  contribution 
of  such  poles  can  be  shown  to  give  an  exponentially  decaying  solution 
in  the  distance  R or  z and  hence  they  contribute  significantly 
to  the  near-field  pressure,  but  their  contribution  to  the  far- 
field  pressure  can  be  neglected. 
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2.4.5  Evaluation  of  Non-Dimensional  Integrals  Contained  in  the 
Far-Field  Pressure  Solution.  To  evaluate  the  non-dimensional 
integrals  (j=l,2,3),  one  can  extend  the  independent  variable  n 
to  the  complex  plane  and  use  the  method  of  Cauchy  contour  integration. 
The  choice  of  the  branch  cuts,  contour  integration  paths,  etc.,  are 
enumerated  below. 

2.4.5a  Choice  of  Branch  Cuts.  The  existence  of  brancli  cuts 
requires  that  the  integrand  be  made  a single-valued  function  before 
Cauchy's  theorem  is  applied  [29].  This  requirement  can  be  fulfilled 


/ 2 2 

by  introducing  cuts  in  the  complex  plane.  Letting  £ = v 1 - q - ■ 


where  n is  a complex  variable,  then  the  sign  of  the  radical 

5 


A - 2 2 


q - n will  be  determined  by  the  radiation  condition 
Im  £ > 0 which  guarantees  the  convergence  of  the  integral  along  the 
entire  complex  path.  The  branch  points  are  defined  by  £ = 0,  which 
yields  the  fixed  values 


/ 2 

= + /I  - q on  the  real  axis.  To  allow 


an  effective  application  of  the  methods  of  operational  analysis, 

)/l-  q2 


the  complex  value  of  (1  + ie/2)i 


has  to  be  used  instead  of 


/l-  2 


q in  order  to  avoid  difficulties  with  the  path  of  integration. 
In  another  word,  at  least  an  infinitesimal  amount  of  material  damping 
should  be  assumed  [27],  Therefore,  branch  points  will  be  located  at 


ie  / 2 

nb  * + (1  + y- ) /I  - q , where  e is  the  elastic  loss  factor. 

After  the  branch  cuts  have  been  determined,  the  undamped  response 
can  then  be  obtained  by  setting  e ■*  0. 

For  the  determination  of  Sommerfeld  branch  cuts,  one  can  start 


with  the  equation 
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r2  L it.  2,  2 

E,  * (1  + 2~ ) (1  - q ) - n 

£ = + i£2  « n = + in2  > 

then 

2 2 

(51  + iC2)  - (1  + f1)  (l  - q2)  - (hj  + in2)2 

2 

= (l  - |-)(i  - q2)  - (ri:2  - n22)  + 2i{|(l  - q2)  - n^}  . 


To  satisfy  the  condition,  Im  £ = 0,  two  conditions  result;  namely. 


£ ,,  2. 

nin2  = 2(1  ‘ q } 


(2.75a) 


and 


(n12  - n22)  < (1  - |-)(1  - q2)  . (2.75b) 

These  two  equations  restrict  the  choice  of  a branch  cut  as  shown  in 
Figure  2.4. 


At  this  point,  one  must  consider  the  changes  of  the  branch 
cuts  necessitated  by  the  limit  e -*•  0 and  choose  the  contour  of 
integration  in  the  upper-half  of  the  complex  r)  plane  either  for  the 
top  Riemann  sheet  (T)  or  bottom  Riemann  sheet  (B)  . The  sign  of  Re£, 
in  the  first  and  second  quadrants  of  the  top  and  bottom  Riemann 
sheets,  can  be  obtained  by  the  following  procedure. 

Let 

5 - /l  - q2  - ,,2  . (/Tv  + )1;J  (/TV  . |'2 

1/7  l°~  + e+  - 71  \ 

= (r_  r+)  exp{  il 2 1}  = yo  exp  { i<J> } , (2.76) 


where  r 


(r  r^1^"  and  <J> 


0_  + 0+  + tt 


are  the  amplitude  and 


phase  angle  of  £,  respectively;  (r_,  9 ) and  (1+»  0+)  the  amplitude 


and  phase  angle  pair  for  vectors  (n 


q ) and  (r\ 


A - q2), 


2 


Figure  2.4  Soramerfeld  Branch  Cuts 
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respectively.  Thus,  the  multivalued  angle  function  0 can  be  defined 
as 

e_  + + 7T 

top  Riemann  sheet:  0 = 

0 + 0+  - TT 

bottom  Riemann  sheet:  $ 

where  the  positive  direction  of  0_  and  0 is  measured  from  the  positive 
real  axis  of  the  complex  plane  p in  the  counterclockwise  sense  as 
shown  in  Figure  2.5.  The  variations  of  the  phase  angle  0 are 
tabulated  below: 


Table  2.1 


Phase  of  Complex  Variable  p 


K2 

AF 

FOD 

DB 


Top 

Riemann  Sheet 


TT 

-r  <_<()<_  tr 

0 1 0 if 

x 71 

* - J 


4>  = 0 

. 71 

0 - J 


Bottom 

Riemann  Sheet 


TT 

2 1 <l>  1 0 

TT 

71  1 ♦ i - 2 


- TT 

TT 

2 


The  above  table  can  be  briefly  described  by 


- nTT  ^ _<  TT  - n7T  , n = 0,  1 , 

where  n = 0 and  n = 1 refer  to  the  top  and  bottom  Riemann  sheets. 


respectively . 
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Examining  the  integrand  in  I ° (j=l,2,3)  together  with  the 
radiation  condition  1m  C > 0,  which  arises  from  the  physical 
requirement  contained  in  the  original  pressure  integral-representation, 
it  can  be  shown  that  the  integrand  of  I ° will  converge  on  the  top 
Riemann  sheet  along  the  entire  closed  path  in  the  upper  complex 
half -plane . 


2.4  5b  Evaluation  of  Branch  Cut  integrals  B.  In  what  follows, 
the  procedures  followed  for  the  evaluation  of  branch  cut  integrals 
will  be  shown  for  the  case  of  the  integral  1^°  only.  The  branch 
cut  integrals  of  1^°  and  I^°  require  similar  procedures  with  only 
slight  modifications. 

The  branch  cut  integral  B of  I9°  may  be  written  as 


B 


+ 


ik  dn  j — t — y 
e ° tq-q^-n  dn 

iaib(n2+q2)  - l)  /l-q2-n2 


+ 


The  signs  of  the  radical  in  the  integrand  can  be  properly 
chosen  by  the  evaluation  of  phase  angle  $ at  the  corresponding 
location  on  the  top  Riemann  sheet: 


r ./  2 2 r 2 “2,  / 2 2 

\1  :/l  - q - h = - / ] 1 - q - ri|=-/l-q  + r 

r / 2 2 /'  2 2|  / 2 2 

r2  :/l  - q - n = rll  -q  — n ! = v 1 — q +r 

r3  \A  - q2  - n2  = - /\l  - q2  - n2|  = - /l  - q2  - r2 

r4  :/l  - q2  - n2  * /\l  - q2  - n2|  » A - q2  - r2  , (2.77) 


where  r is  a positive  real  variable  representing  the  radius  of  the 


corresponding  vector. 


48 


where  the  integral  of  (F^  + F^)  is  first  divided  into  two  subintegrals 
(from  0 to  1 and  from  1 to  ®),  and  then  the  subintegral  from  1 to  00 
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is  again  retransformed  to  an  integral  from  0 to  1 for  numerical 
evaluation . 

2.4.5c  Evaluation  of  the  Pole  Contributions.  The  positions 
of  the  poles  are  determined  by  the  zeros  of  the  denominator  of  the 
integrand  in  the  integral  1 ° (j=l,2,3),  i.e., 


i{b(n2  + q2)  - 1}  A 


q2  - n2  +1=0  . (2.79a) 


If  one  lets 


A 2 n2 

/I  - q - n 


q - n = is  , 


Equation  (2.79a)  becomes 


s5  + 2s3  + s(l  - ~)  - y = 0 , 

* 


where 


(2.79b) 


* ■ ^ . 

p UJ  o 

o 

Once  the  roots  of  Equation  (2.79b)  are  found  for  different 
frequencies,  the  variable  will  be  obtained  in  terms  of  the  variable 


s by  the  following  relationship: 


/ 2 A 2 

/I  - q + s 


(2.80) 


From  Equations  (2.80)  and  (2  79a),  it  becomes  obvious  that 
there  exists  ten  poles  in  the  integrands  of  I ° (j-1,2,3).  However, 
only  those  roots  on  the  top  Riemann  sheet  are  needed  for  the  use  in 
the  Residue  theorem.  The  roots  that  fall  in  the  top  Riemann  sheet 
can  be  found  from  the  mappings  between  variable  s and  imposed 
by  the  radiation  condition  Im  £ > 0.  Since 


A - q2  - h2  = is  , 
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or 

^1  + i^2  * i(®R  + isI)  = ~ si  + isR  » 
the  condition  Iin  C > 0 is  seen  to  correspond  to  the  condition  Re  s > 0 . 

Numerical  computation  shows  that  the  typical  five  roots  of  Equation 

(2.79b)  are 


1 positive  real  root:  (sD1) 


4 complex  roots:  (sR2  + isI2)  » (sR2  “ 

(-sR3  + isI3)  , (-sR3  - isI3)  , 

where  sD,  (j=l,2,3)  and  s.,,  (j=2,3)  are  positive  real  numbers. 

KJ  ij 

Thus,  the  desired  roots  in  t|  can  be  found  first  by  emp] eying 
the  condition  Re  s > 0 to  pick  out  the  three  roots  in  s which  consist 
of  one  positive  real  root  (sR^)  and  two  complex  roots  (+sr2  + isl2), 
(+s^2  ~:*‘sx2^  anc*  t^en  substituting  these  into  Equation  (2.80)  to 
obtain  six  roots  in  r|  on  the  top  Riemann  sheet.  Out  of  these  six 
roots,  only  the  three  roots  which  include  one  positive  real  root  and 
the  two  complex  roots  lying  in  the  upper  half  plane  of  the  complex 
variable  n are  the  desired  ones.  The  reason  for  excluding  the  negative 
real  pole  is  that  if  a slight  amount  of  damping  is  assumed,  this 
negative  real  root  will  be  shifted  off  the  real  axis  into  the  lower 
half  complex  plane  1 and  outside  the  domain  of  the  closed  contour 
integration. 

2.4. 5d  Computation  of  Residues  R . In  general,  the  residue 

3 

of  A(n)/B(n)  at  a simple  pole  n = r)^  can  be  found  by  the  formula: 

d . Mn) 
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If  the  expressions  A(H)  and  B(ri/  contain  multi-valued  functions,  such 


as  a square  root,  the  signs  of  the  radical  have  to  be  carefully 
evaluated.  Thus,  if  one  takes  the  integral  1,,J,  lor  example,  t 

residues  of  this  integral  are: 

ik  dn  . — - — - 
o / 2 2 

R e /1-q  -n 


he 


|^{iaVl-q2-n2(b(n2+q2)  -1]  + 1} 


n = ir 


ik  dn  . — — 

o / 2 2 

e /1-q  -n 


ia{4br)(n2+q2)*/l-q2-n2  - n[b(n2+q2)  -lj(l-q2-n2) 


-1/2 


n = n 


ik  dn 

o j 2 2, 

: -.(.l.-.q  -Hj.  .) 


ianj{b(nj2+q2)(-5nj2-5q2+4)  + 1} 


, j-1,2,3 


It  is  noted  that  the  numerator  and  denominator  of  the  residue 
expression  is  simplified  by  multiplying  both  by  the  same  branch  of 


A - 


- rr 


, resulting  in  a final  residue  expression  which  is  free 


from  the  multivalued  radicals. 


2 4.5e  Final  Solutions  of  the  Integrals  1^  (j=l,2,3).  Using 

the  methods  of  evaluating  branch  cut  integrals  and  residues  developed 
in  previous  sections,  the  integrals  I ° (j=l,2,3)  can  be  evaluated 
according  to  the  theory  of  contour  integration  as  follows: 


- - 

i/J  ♦/  ♦/! 

( 

. [Integrand]  + 2ru/ 

f>  \ 

) Residue) 

J J j \ 

r r,  r„  r 

1 1 

r1 

1 4 '3  4 

J 

or 


'30 

‘i°  ■ -2/ 


A 2-  2 


q +r  d 


2 2 1 +2/  ' 2 2 

o 7 7 7 7 J > ~>)  7 7 

0 l+a  fb(q  -r  ) -1}  (l-q^+r  ) 0 l+a" fb (q~+i “ ) -1}  (l-q^-r  ) 


Arq  / 2 2 

/ 1-q  -i  dr 
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£ d-q2-nn2) 

2TT!  £ 

1 iann(b  (nn2+q2)  (-5Hn2-5q2+4)-l} 


oo 

■-21/ 


* -k  dr 


ur  / ~ s- 

° / 1-q  + r"  d 


v'  i 

— +2  y 


-q  lk  dr  / ~ ~ 

0/22, 
e /1-q  -r  dr 


I 2 2 2 ^ ^ 22  J 2-  22^  ^ 22 

0 l+a  {b(q  -r  ) -1}  (1-q  +r  ) 0 1+a  {b(q  +r^)  -1}  (1-q  -r  ) 


+ 2rri  V* 
r^l 


ik  dn  9 2 

e 0 W-nn) 


iann{b(nn2^q2)(-5nn2-5q2+4)-i} 


A-a2 


® -2k  dr  / — 5 5-  /1-q  2ik  dr  / 5 9 

r 0/2  2 r o/22, 

I e v 1-q  +r  dr  / e /1-q  -t  dr 

" / 22  +/  22 

0 l+a  {b(q  -r  ) -1}  (1-q  +r  ) 0 l+a^{b (q  +r  ) -1}  (1-q  -r) 


2ik  dnrt  oo 

£ e ° n (l-,2-nn2) 

n-1  lar\„(b(ntl2+q2)<-5nn2-5q2t4)-ll 


2.4.6  Approxlmace  Expressions  for  the  Integrals  1^°  (j=l,2,3). 
Although  a digital  computer  can  be  used  for  the  numerical  computation 
of  the  integrals  I ° (J*l,2,3),  which  consist  of  branch  cut  integrals 
and  residue  terms,  the  resulting  numerical  values,  when  substituted 
back  into  the  radiated  pressure  in  Equation  (2.74),  do  not  completely 
allow  one  to  interpret  the  trend  in  the  variation  of  radiated  far- 
field  pressure.  Thus , in  order  to  gain  more  insight  into  the  behavior 
of  the  integrals,  it  seems  plausible  to  obtain  the  approximate 
expression  for  the  integrals  I^°.  These  asymptotic  expressions  allow 
a better  understanding  of  the  fundamental  characteristics  which  govern 


the  physical  phenomena  of  the  radiated  far-tleld  pressure  distribution. 
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2. A. 6a  Approximate  Expressions  of  the  Integral  I^°.  To 
begin  with,  one  can  express  I 0 m terms  of  branch  cot  integrals 
and  residue  contribution  as  developed  in  previous  sections,  i,e  , 

T O _ f /i-q“-n2  dn 

1 = / H 5 = A " 1(B+C> 

J / 2 2 9 2 L 

-°°  ia/l-q  -n  b(n  +q  ) -l}+  1 
1 , 2 2 
♦ y 1-q  ~nn 

^ 2 2 2 2 > 
n=l  iann{b(nn  +q  1 (-5n n -5q  +A)+1 } 

where 

, - ■“  9 *•  ■ ’ • .0a.  b .(M2  • «2  , 

O o o 


AAA 


22 
i/l-q  -r 


J 2 2 2^2  22 

0 1+a  {b(r  +q  ) -1}  (i-q2-r2) 


2/l-q2+r2  dr 


J o 22^^  29 

0 l+a  {b(q  -r  ) -1}  (1-q  +r  ) 


, / 2 

1 / /1-q 


2r  t'/l+(l-q2)r2  dr 

r10+a2{b<q2r2-l)2-r4Ml+(l-q2)r2i 


2 ( 1-q  ) 


r?/l+r2 


0 r10+a2(l-q2){b[q2r2-(l-q2)  f-rVd+r2) 


(1)  Low  frequency  approximation  (ft  <<  1) : 

y/C- q2  


A = 2 


/ 2 2 , ,,  2,  Tt 

/1-q  -r  dr  = (1-q  ) j , 


/i-„2 


' 2j  >4-q2+r2  dr  “ (l-q~){/2  + log^Q  + /2)}  , 
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and 


C = (l-q“) 


r V 1+r  dr 


J ^ 2 2 2 999^ 

0 r +a  (1-q  )(b[q >(lV)l  -i  ) (1+r) 


(1-q  ) ft 


“ <W H 


0 


5 2U 

r +a  b 


r3  dr 

b2 (1-q2 ) 5 I 


(1-q  )ft 


1/5  2 

U (l-q2)ft  , 

• aV»/f  Tf  * | 

r U 91/5 


U 

, , 4 9 14 

/,  2.  r r x 

(lq)4u“2  3 

9 U I4u 


r +p 
. 1/5 


(1-q  )ft 


6r  Hr 


11 


+ -"J  +5 

.1/5 


^ "2/5)02“6/5)  = 1 . 135 (a  _2/5)(ft)'6/5 

L I I L O O 


[Note  that  u = a2b2(l-q2)  ] 

d i ~ -1/5  n~3/5 

Poles:  a0  (2 

n,  a°-1/5  sT3/5  e2"1'5 

2 o 

„ _ -1/5  0-3/5  3^1/5 

21  a ft  e ' 

J o 


Summing  up  the  various  contributions  evaluated  above,  one  obtains: 

22  (1-q2)  | - i{(l-q*')[/2  + Xoge(l+/5")]  + 1.135a  -2/5  fi'6/5}+r 

n^l 

(2)  High  frequency  approximation  (ft  > 1): 

(i)  When  q2  < ^ « a: 


Res 


55 


A = 2 (1-q4) 


f /l-t  dr _____ 

1 . 2,  2 2,  ,,  2.  . , . 2,  2 2,' 2 2,2 
0 1+a  b (1-q  )(l-r  )[{(l-q  )r  +q  } -a  ] 


- a-O  [V3? 2 

•'o  ' l+a2b2(l-q2)(l-r)t{(l-q2)r+q2}  -C2]2 


a f dr 

2 Vrt-  2 I , , , 2 2 

°‘q  -»l  + aV(l-o)[{(lV)r}  -a2] 

. Aiaip^i °v  ■ ■ 

TO- -q2  4 (1-q2)  ( 1-0 ) io3/a  1+T2 


TO  Ttfi  . 1 , 

= — ■ , (Note  that  T = ) 

0 2 0 / /f)  2 2a  /Iro 

2a  /O-q  2a  /l/S2-q  o 

o o 


- (1-q  ) 


2 (1-q2)  f 

Jo 

CO 

JV1? 


/l+£2 


j /1+r  dr 

I 2 2 

f 29  9 9 999  9 

0 1+a  b (1-q  )(l+r  )I{(l-q  )r  -q  } -o  ] 


22  2 2 222 

1+a  b (1-q4) (1+r) [(1-q  ) r -o  ] 


llfL.  f dr 

2 V 2^t  / , , , 2 , _ 2 

M 1+a  b (1+a)  [(1-q  ) r -a  J 


- / 2 

2a  TO+q 
o n 


Poles:  n.  - /l-q2 


/ 2 

+ l/o+q 


/ 2 

4a  J 1+0  /a+q 


3 /— J 

4a  / 1+0  /o+q 
o 

3 2 

T+  ..  1TO 

2 Res . - 

n-1  a /+q2 

o 


+ it'O+q 


. / . 2 

i/O+q 


/ . 2 
- l/o+q 
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Summing  the  various  contributions  obtained  above,  I 0 becomes: 


2 2 
T ° . , iTTO 

^ + _ 


2a  «/o-q2  2a  A+q2 


(ii)  When  q^  = 1/ft  = o , 


A = 2 ( 1-q  ) 


A-t2  dr 


J 22  2 2 1 7 o 2 o 2 

0 l+a  b (1-q  ; (1-r  ) [{ (l-q2)r  +q2}  -O2] 


- 2(l-q  ) 


0 l+4a2b2q4(l-q2)  r4 


» (1-q  ) 


2,  2 4,,  2/  4 

l+4a  b q (1-q  ) r 


(l-l/fl)*'*  <1 


4abo (l-o) ' 


TTO(l-O)' 


4abo (l-o) 


no (l-o) 


4abo (l-o) 


I , 2 2 

Note  that  there  are  two  critical  points  |r  -13—  , 3_~q 

1 , 2 ,lZ  . 2 

' 1-q  1-q 

existing  in  the  integral  D. 


. 3/2 

£ Res.  , 122 

n = l /2  a 

o 


Summing  up,  the  integral  1^°  becomes: 


0 , nqu-q)1/4  /, 

1 V'2  ' 


4 abo (l-o ) 


}/7  I 2..1'2  V 4abo  ( l-o ) ^ 2 


(iii)  When  q2  * 1/ft  = a: 
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2 ( 1-q 

2 2 2 4 

1+a  b (1-q  )(q 


J. 

1 [ 

-oV  l 


A-t‘ 


dr 


no 


o 2<  4 „2\2 

2a  (q  -O  ) 

o 


I 2 2 

■>'=-!- 


o / 2 / 2 

✓o+q  /q  -0 


3 . 2 

_ 1TO 

Res.  - — 

a A+q 
o 


Therefore , 


o .. 


TTO 


in 


1 2 4 2 2 2a 

2a0  (q  -a  > 


Jo\  / 2 fl  / 

\ /q  -a  /q  +0  / 

2.4.6b  Approximate  Expressions  toi  the  Integral  1 
1 2°  can  be  expressed  in  the  following  form: 


°°  ik  dn  / — ~ — o 

° /l-q2-r,2  dr 


/ 2 n 2 2 “ 2 
-“>iab/l-q  -ri  I ( n +q  ) -o  ] + 1 


= A - lD 


3 

I 

n*  1 


lk  dn  „ „ 

o n 2 2. 

e (1-q  -nn  ) 


iari^b(n|2+q2)(-5r^2-5q2+4)  + it 


(1)  Low  frequency  approximation  (ft  ' 1) : 


(i)  When  1/2  < k d < 1: 
o 


/l-q‘ 

J 


lk  dr 


2e 


cl  -r  dr 


2 2 
?2  22  22  2 

1+a  b (1-q  -r  )t  (r  +q  ) -0  ) 


Integral 
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/ 2 

ik  dr  j f-q 

/l-q'-r2  dr  2 f O 

0 7) 


a-q 

2 1 e 


i 


He  dr) 

— + lk  dr' 

2 o 


/ 2 2 . 

l-q  -r  dr 


. (l-q‘)7T  2 2. £ tt  2kndl  7 

2 (k  d)  > 7Z  + “4-  d-q  > 


3/2 


16  3 


D a 2 (1-q  ) 


-30 

/ 


, ,/  2 

-k  d/l-q  r / 7 


/l  2 


1+r  dr 


2 ( 1-q2) ^ / e ° 


0 l+a2b2(l-q2)(H-r2)[{(l-q2)r2-q2}2-  g2) 

1 -k  dr/l-q-  i 2l  r -kdr/l-q2 


/l-q2 


-k  d/l-q 


/l-q2 I 


21  r -k 

1+y-ldr  + f e 


-k  d 


A-o2 


k d 
o 


1-e 


Poles:  r),  = a ~1/5  ft  3/5 
1 o 


n - a “1/5  o“3'5  27T1/5 

Ho  - a u e 

^ o 

_ ~ -1/5  0—  3/5  3m/ 5 

n,  - a ft  e 

J O 


r dr] 


Thus,  the  integral  has  the  following  approximation: 


T o (l-q2)n  ^V)2(l-q2)2  2kodi(l-q2) 

2 2 ' 16  + 3 


3/2 


- 1 


/l-< 


2 2-e 


-k  d/l-q 


2 1-e 


-k  d/C/) 


k d 
o 


/-i  k V 

2 1-q  > o 


3 

+ £ Hes  . 
n-1 


lil)  When  k d > 1: 

o 


/ ? 

✓l-q2 


lk  dr 


° /,  2 2 . 

/1-q  -r  dr 
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2 ( 1-q  2 ) / e 


lk  dr21-q2 


/i  r.  ui»  x-u  | 

u /— v ' 


1-  y-l  dr 


= 2(l-q  ) 


ik  dr'/l-q  ~ lk  dr/l-q  lk  dr/l-q 

o 2 o n o 

e re  e 

+ 

Ltk  d/l-q2  2 ik  d/l-q2 

o o 


iX-q2) 


ik  dr*4-q2  - 1/ 

I ° ( 


/ J 


(ik  d/l-q 
o ^ 


0.r,  2,2,,  2 ik  d^-q2  . , 

2i’k  d ( 1— q )-l)  ;e  ° / / — ^ lk  d ( 1— q ) 

+ /-k  dA-q  -l-  — — =- 

o 2 


2,2,,  2, 
: c 
o 


V/-„2 


k 3dM-q2 

o 


D 


k d 
o 


(obtained  from  case  (i) ) . 


Poles:  r|.-  3/;i  ft  >>l  3 (dominate) 

r),~  a°-1/5  ft~3/5  e2lTl/5 

2 o 

-1/5  n-3/5  3iri/5 
H_-  a ft  e 

3 o 


It  is  not  difficult  to  understand  that  the  branch  cut 
integrals  "A"  and  "D"  may  be  negligible  compared  with  the  residue 
contribution  due  to  the  real  pole;  thus. 


I 


Res . 


(2)  High  frequency  approximation  (ft  > 1): 


(i)  When  q < 1 /ft  - o: 


A ~ J 


/.  2 

vl-q 


2e 


lk  dr 
o 


A- 2 2 


q -r  dr 


. 2.  2..  2 2,  r , 2 2 2 2 , 2 

1+a  b ( 1-q  -r  )((r  +q  ) -O  ] 


2 ( 1-q  ) 


1 

L 


t/l-q2 


ik  dt/l-q 
o 1 


A-rl 


dr 


no 


0 l-»-a2b  2 ( 1-q  1- 1 2 ) [ { ( 1-q  2 ) r 2+q  2 1 -O2] 

2 ik  d/o-q2 

r\  * 


/ 2 
2a  /0-q 
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-k  dr/l-q 
o 


D * 2 (1-q" ) 


/l+r2  dr 


) y 2 ? 9, 

0 l+ab(l-q)(l+r)[{(l-q)r-q)  -0  ] 


.2  -k  d/0+q 

TO  O 

e 

/ 2 
2a  /o+q 
o 


The  expression  for  "D"  is  expected  to  be  a good  approximation  under  the 
conditions  of  high  frequency  (ft  >>  1)  and  not  very  large  k d.  In 


another  word,  the  approximate  result  is  good  for  when 


-k  dr/-q2 

e is  a slowly  varying  function. 

3 . 2 -k  d/+q2 

£ Res.  * e ° 

n=1  2a  /o+q2 

o 

Thus,  the  approximate  expression  for  I^0  can  be  obtained  as  follows: 


O TO 


2 ik  dJo-q- 


1 /o+Q  2 


/ 2 

/a-q 


2 d/a+q^ 

i7T0  O M 

+ e 


2a  /a-q 
o 


/ 2 
/O+q 


2a  /o+q 
o n 


(ii)  When  q « 1 /ft  = o: 

(a)  k d/l-q2  < 1: 
o 


2 ik  dr/.-q'- 


1 n £_\  ° s 

2 r 2 2 e 

A * 2 (1-q  / / - 

J 224  2.4 

0 l+4a  b q (1-q  ) r 


Since 


ik  dr/l-q2 
e 0 » 1 - 


k d (1-q  )r  /— 2 

-2 + ik  dr/l-q'  , 


Similarly,  since 


A « iqlAzg.il7.4  A i \ 

2ao1/2  j 4abo(l-o) 2/2 


TO3k  2d2 
o 


„ 3/2,.  ,1/4/  4a 

8aQ  (l-o)  I o 


ink  do 
o 


-k  drZ-q2  y — j 

e ° - 1-k  dr/l-q  , 

o 
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D = 2 (1-q  ) 


'/ 


-k  dr/l-q2 
o M 


/l+r2  dr 


22  2 2 222^2^ 

0 l+aV(lV)(l+r  )UU-q  )r  -q^}  -a2] 


= (1-q  ) 


(i-k  dr/l-q  ) dr 
o 2 

, ...  2.2  4.,  2,3  4 

' „ l+4a  b q (1-q  ) r 


/ 


,1/4  nk  da2  3/2  -k0d,/,2o 

TO  (l-o ) ^ u TO  e 

2a  1/2  4ao  2/2  a 

o c 


3 3/2  -k  d/2o 

r-  « _ lTO  o 

) Res.  = e 

n=l  /2  a 


Thus,  the  integral  I,  can  be  represented  by: 

.14/ 


O TO ( 1-0 ) ' 

2 “ 1/2 


1 - 


3.  2,2 
TO  k d 
o 


2a 


4abo ( 1-0 ) 3/ 2 f 8a  3/2(l-o)1,4 

o 


/,o(i-o)1/4  . 

1 1 2a  1/2  2a  " 2 JTT  ^ 

o o o 


(b)  k d/l-q  > 1: 
o 


r/l-o2 


A = 2(l-q  ) 


. 2 ik  dr/l-q 

l/i  r , o M 
(1  — e 


r 11  " T 

/ — i 


dr 


2 2 4 2 4 

0 l+4aS  q (1-q  ) r 


(1-q  ) 


00 

/ 


/l-c2' 


cos(k  dr/l-q  ) dr 
o 

...  2.2  4..  2, 3 4 

1+Aa  b q (1-q  ) r 


+ 2 ( 


1-q2)  / 


1 / 2 

sin(k  dr/l-q  ) dr 
o 


„ 1^,  2,  2 4,.  2,  4 

0 l+4a  b q (1-q  ) r 


-sin 


TOk  d/4 

o 


TO(l-O) 
/T~  a 


1/4 

lTT 


/2-  a 1/2(l-o)1/4 


TOk  d/4 

.1  0 it 

sin| cos  T77  + 71 

l/2rt  „\l/4  4i 


/2a  ^ClTO) 
o 


/ 


- 


62 


0 

2, 

g k d , i 

+ ° tan"1 

^2a(l-o)’/2 

1 2 (k  dg)2  [ 

1 

l 0 J 

and 

2(l-o) 


1/2 


k d 
o 


3 ,3/2  -k  d/2<7 

£ Res.  = e ° 

n=l  /2a 


-sin 


Tick  d/4 
o 


1/4 

T o . rrg(l-g) 

2 

o 
2 

Ok  d 


/2a  1/2(1 -o)1/4 


sin 


co  = 


TOk  d/4 
o 


/2a  1/2 (l-o) l' 4 
o 


4 


+ t -iV.d-g)122!  2i (l-o) 122  ^ i„o3/2  -K0^ 

' 2(kodo)2  ( kod  /2a 


(iii)  When  q > l/fl  = o: 
(a)  k d/l-q2  < 1: 


A = 


■lllzl  ). 


u 2,2. 1 2,.  4 2. 

1+a  b ( 1-q  )(q  -a  ) 0 


-L 

I 


ik  dr/l-q2 
i /l-r  dr 


TTO 


i/l-q2(g6) 


+ i( 


1 2k^d/l-q  (g  ) 


D = 


-2,4  2 2 'j.  2.  4 2 2 

2a^  (q  -O  ) ^3ao  (q  -g  ) 

/ 2 

_ i _ 2 -k  dv'O+q 

tt  y a o M 

2a 


° t/a‘,  2 


2 -k  d/q 
g o M 

e 


/ 2 

/q  -g 


Jl 

/q  -C 


/g+q  /q  -g 

Ihe  expression  for  ' D"  is  expected  to  be  a good  approximation  for 
high  frequency  (Jl  >>  1)  and  not  very  large  values  of  k d, 


I Res. 

n*l 


ino2  -R0d^+q2 
e 


/?  2 
a /g+q 
o 
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Thus,  an  approximation  for  I ° can  be  written  as: 


I 


7T0 


, 2,  4 2 2 

2aQ  (q  -o  ) 


2k 

+ i<— £ 


d06(l-q2) 


l/; 


T 2,  4 2 ‘ 

3aQ  (q  -o  ) 


in 

2a 


-kod/a+q2  o2  -kod/q2-a 


/,  2 
/o+q 


/ 2 

/q  -a 


(b)  k d/l-q 
o M 


/l-q2  > 1: 


A = 


■ ( l-q2 ) 


l+a2b2(l-q2)(qA-oV  "0 


X 

/ 


lk"dr/l"q2  £7 


dr 


2(lV) 


2 2^-1^  a 2 ~ 

a (l-q")  (q  -O)  (k  d)J 
o o 


i { k 2d2(l-q2)-l}  + 
o 


/ 2 

lk  d/l-q 


2 2 2 
/— 2 ik0  d' <W > 

-k  d/l-qZ  - i 5 

1 o 2 


and 


„ _ n / o2  -k0dv4+q2  2 -knd^2-o 

D-— <— e + e 

/ 2 

/q  -a 


° U 2 


/O+q 


llie  expression  for  D is  valid  for  high  frequency  (£2  ■>  1) 


and  not  very  large  values  of  k d 

o 


, 2 -k  d/^j+q" 

£ Res.  = — - ° 


3 

n*l 


2 

a /O+q 
o n 


and 


4o 


2 9 1/2  A o 2 T 

a n (l-q  ) (q  -a2)  (k  d)3  L 

o o 


/ 2 

99  9 I d/l-q 

i^kVd-q2)-!^  e 0 
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/— 2 Ik  2d2(l-q2)l 

- k d/l-q  - i 2 — r 

o i 


ITT 

2a 


„ -k  dv'O+q  . -k  dv' q -a) 

, 2 c 2 o 

0 e 0 e 


/ 2 

/o+q 


/ 2 

/q 


q -0 


2.4.6c  Approximate  Expressions  tor  the  Integral  1^  . The 

solutions  of  the  integral  I^°  can  be  obtained  by  replacing  (k^d) 

by  (2k  d)  in  the  solutions  1„°. 
o d 

2 . 5 Evaluation  of  the  Radiated  Acoustic  Power 

2.5.1  Acoustic  Power  Generated  by  Far-Field  Pressure  (N^) 
The  radiated  acoustic  power  can  be  evaluated  by  integrating  the 
square  of  the  radiated  far-field  pressure  over  a large  spherical 
surface,  as  follows: 


tt/2  2tt 


ff 


* // 

A-n  a. 


pvR  sin  <p  d(p  d8 


=0  6=0 


2P„ 


tt/2 

r 

2tt 

r 

L 

N 

O 

7 

tt/2 

j 

4-0 

J 

0=0 

!P(R,  <S>,  9)1  sm  $>  dT  d9 


| p(R,  <p,  0)|  sin  <J>  d<J>  dQ  , 


where 


p(R,  4>,  0) 


ik  R 

-ik  F e ° cos  2ttR 
o o 

ik  m i 2 . 

, O J,  u)  J 4 A 

1 cos <p|  1 r-  sin  <p 

p I 2 

o l u) 

o 


C (0,  <f>) 
ms 


and 


65 


C (6,4)) 
ms 


1+U, 


mi  1/  W1  > 

~\l  2 5 — I Cos(k  6 COS  0 COS  4>)  - I ° 

Poh/ [a:,  -w  / 1 ° l 


- I. 


1 - 11,1, 


o J M 


jPoh  IpJ,  -u 


2 II  2“  U'L  " 


/ ml) 

Iph2 

2 j M 

J1  ,2 
\pQh 


mi  w ^ 

,2  2 2 

P h w.  -u) 


+ u3uy  *ln‘  * sl"‘  4i°  - “il-TlII-lH 

Ip  1)  l\  U).  -CD 
o 1 


v ♦ w - v 


For  the  complete  definition  of  "C^,"  one  tan  refer  to  Equation 
(2.74).  Since 
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the  radiated  acoustic  power  can  be  written  in  the  following  integral 
form: 
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(2.81) 


2-5.2  Acoustic  Power  Generated  at  the  Surface 
of  the  Plate  (N  ) . The  Fourier  transform  of  the  velocity 
distribution  over  the  surface  of  the  infinite  plate  under  discussion 
can  easily  be  obtained  by  substituting  Equation  (2.57)  into  Equation 
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The  bracketed  inner  integrals  have  the  properties  of  a two-dimensional 

Dirac's  delta  function  6(k  - k ')  6 ( k - k ')  which  can  be  derived 
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from  the  following  Fourier  transform  pair: 
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Thus,  the  expression  for  the  radiated  power  can  also  be 


written  as: 
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where  the  relation  p(k  , k , z=0) 
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substituted  to  arrive  at  the  above  result. 

The  expression  for  the  radiated  power,  given  in  Equation 
(2.82)  can  be  simplified  under  the  transformations  of  independent 
coordinates 
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Hence,  Equation  (2.82)  becomes 
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The  limit  of  integration  on  k is  extended  only  to  k = k , because 

r 1 r o 

for  k^  > kQ,  the  real  part  of  the  integrand  vanishes. 

One  can  further  simplify  the  expression  for  acoustic  power 
k 

£ 

by  substituting  sin  <J>  = , dk  = k^  cos  ' d<$>'  and  noting  that  the 

o 

response  of  the  entire  structure  is  symmetrical  about  the  x and  y 
axes,  i.e.,  the  integrand  will  remain  unchanged  when  (0')  is  replaced 
by  (-0 ' ) or  by  (tt— 0 ' ) , 

Thus,  Equation  (2.82)  becomes 
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It  should  be  noted  chat  the  expression  for  the  radiated 
acoustic  power  given  in  Equation  (2.83)  is  exactly  equal  to  that 
given  in  Equation  (2.81).  This  equivalence  pi  ides  an  analytic 
verification  for  the  asymptotic  solutions  obtained  for  the  far-field 


acoustic  pressure. 


CHAPTER  III 


NUMERICAL  RESULTS 


3 • 1 Numerical  Results  for  a Beam-Reinforced  Plate 

3,3  1 Input  Data.  For  an  aluminum  infinite  plate  m 
water,  numerical  results  were  obtained  for  the  following  physical 


constants : 


E = 10  x 10  psi  ; V = 0.33, 


p (density  of  plate)  = 2 56  x 10  4 lb-sec~/in4 , 

Po  (density  of  fluid  medium)  = 0.948  x 104  lb-sec2/ln4 


c = 5.64  x 10  in/sec  ; R = 1 yd.  = 36  in. 


3,l-2  Far-Field  Pressure  Spectrum  and  Directivity  Functi.  n. 
The  solution  for  the  far-field  pressure  for  the  point  excited  beam- 
reinforced  plate,  as  shown  in  Figure  2.1, can  be  derived  directly  from 
Equation  (2.74)  by  letting  m,  0 and  -*•  0 (or  to  0)  The 
reduced  expression  becomes: 
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~2  ; m = p^h  (mass  of  plate  per  unit  surface  area); 
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The  expression  within  the  bracket  in  Equation  (3.1)  represents 

the  radiated  far-field  acoustic  pressure  of  an  infinite  isotropic 

plate  due  to  a point  force  excitation.  The  nond  linens  lonal  coetficient 

C stands  for  a pressure  amplitude  modification  factor  when  the 
ms 

plate  is  reinforced  by  a line  mass,  a stiffness  member  or  bean-like 
discontinuity.  The  excitation  is  delivered  by  a point  force  located 
at  the  discontinuity.  This  coefficient  is  a function  oi  the  mass  and/or 
the  stiffness  of  the  attached  reinforcement,  the  excitation 
frequency,  the  plate  material  properties,  the  sound  velocity  as  well 
as  the  location  of  an  observation  point.  The  dimensionless  variables 
and  both  depend  on  the  excltati>  a frequency  and  the  material 
constants  of  fluid  and  plate.  It  is  worth  mentioning  that  p is  a 
function  of  the  excitation  frequency  to  the  second  power  and  is  related 
to  the  mass  effect  of  the  attached  reinforcement,  while  ^ is  a 
function  of  excitation  frequency  to  the  fourth  power  and  is  associated 


12 


with  the  stiffness  of  the  attached  reinforcement.  The 
dimensionless  variable  when  multiplied  by  the  fourth  power  of  the 
thickness  of  the  plate  represents  the  effective  moment  of  inertia  of 
the  attachment.  The  nondimensional  integral  1 ° indicates  the  strong 
coupling  between  a beam-reinforced  plate  and  an  adjacent  fluid.  The 

Q 

integral  1^  strongly  depends  on  the  excitation  frequency  as  well 
as  the  location  of  an  observation  point.  The  coincidence  frequency 
to  is  the  frequency  at  which  the  flexural  wavelength  in  the  plate, 
as  predicted  by  the  classical  plate,  equals  the  acoust  c wavelength. 

It  should  be  noted  that  the  thickness  of  the  plate  h appears 
explicitly  in  the  coincidence  frequency  in  only;  hence,  the 
following  computations  were  made  such  that  h does  not  appear  explicitly 
in  the  expression  for  the  far-field  pressure. 

With  Equation  (3.1)  in  hand,  the  following  numerical 
computations  were  made: 

(l)  Computations  were  performed  for  the  prediction  of  the 
far-field  spectrum  at  zero  azimuthal  angle  (4>  = 0J)  per  1 pound  force 
at  1 yard  in  db  relative  to  1 (Jbar-in  versus  normalized  frequencies 
(£1)  at  five  different  mass  ratios,  i.e.,  y * 0,  1,  2.7,  3,  10. 
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The  pressure  spectra  are  plotted  for  the  cases  of  Y»0,l,10 
shown  in  Figure  3.1.  The  curves  for  y<*  2.7  and  5.0  are  not  shown. 
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although  they  fall  between  the  curves  y = 1 and y = 10.  It  should  be 
noted  that  the  stiffness  of  the  attachment  does  not  have  any  influence 
on  the  pressure  spectrum  at  this  particular  point  of  observation. 

(ii)  The  directivity  function  for  the  radiated  far-fieid 
pressure  from  a point-excited  beam-reinforced  plate  was  computed 
relative  to  the  far-field  pressure  at  zero-azimuthal  angle  for  a point- 
excited  isotropic  plate  versus  the  azimuthal  angle  <)>  at  different 
rotational  angles  (0=0,  90°  and/or  45")  and  different  excitation 
frequencies . 
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Directivity  computations  were  made  at  the  frequencies  (2  = 0.01, 
0.1,  0.5,  1.0,  1.5*''*2.0,  2.5  for  the  following  three  types  of  the 


reinforcement;  namely, 

(1)  Mass  reinforcement:  y = 0,  1.0,  2.7,  5.0,  10.0 

(2)  Stiffness  reinforcement:  Since  the  stiffness  of  the 

reinforcement  can  be  characterized  by  the  ratio  of  the 
depth  of  the  attached  rectangular  beam  (qnh)  to  that  of 
the  plate  (h)  while  assuming  the  width  of  the  beam  is 
equal  to  the  depth  of  the  plate,  five  different 
values  of  stiffness  are  considered,  i.e., 

q2  - Ol-p-  = ol,  0 . 37  (y^—  * 0.051  , l.of-p-  = 0.89 

' P ' ' P ' * P 
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where  I = bD  = hD  , 
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V = 


3 3 4 

b(q2h)  q2  h 


(Moment  oi  inertia  of  the 


12  12 
attached  beam  about  its  own  neutral  axis). 

(3)  Beam-like  reinforcement:  {y  = 0,  q2  = 0}, 

{y  = 1.0,  q 2 = 0. 37} ,{y  = 2.7 , q2  = 1.0}, 

{y  = 5.0,  q2  = 1.85}  (y  = 10.0,  q.,  = 3.7}. 

Plots  of  the  directivity  function  were  made  for  (2  = 0.5,  1.5, 
and  2.5  only.  The  case  of  SI  = 0.5  has  been  chosen  for  a representative 
excitation  frequency  which  is  below  the  coincidence  frequency,  since 
the  results  for  very  low  frequencies  such  as  S2  = 0.01  and  0.1  do  not, 
in  general,  show  pronounced  variation  for  moderately  sized 
reinforcements.  The  case  of  S2  = 1.5  shows  the  results  slightly 
above  the  coincidence  frequency,  which  not  only  represents  the 
phenomenon  around  the  coincidence  frequency  but  also  shows  more 
definite  and  noticeable  variation  than  the  case  of  ft  = 1.0.  The 
frequency  at  SI  = 2.5  is  representative  of  frequency  of  excitation 
much  higher  than  the  coincidence  frequency. 

It  was  noted  earlier  that  the  pressure  spectrum  computation 
at  $ = 0"  was  not  influenced  by  the  stiffness  of  the  attachment. 
Similarly,  the  pressure  spectrum  at  the  particular  rotational  angle 
♦ a O’  is  unaffected  by  the  stiffness;  hence,  a reinforcing  beam 
acts  just  like  that  of  its  mass  only  for  the  directivity  plot  at 
<t>  « O’. 


There  are  seven  figures  plotted  for  Section  (ii): 

(1)  Below  the  coincidence  frequency  at  il  » 0.5:  Figure 

3.2  was  made  at  two  fixed  rotational  angles  (0  * 0°  and 
90°)  for  various  types  of  reinforcements. 
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Figure  3.2  Directivity  Function  of  Far-Field  Pressure  of  Beam-Reinforced 


(2)  Above  the  coincidence  frequency  at  (2  15:  Three 

figures  were  plotted  according  to  three  types  of 
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reinforcement  Figure  33  was  made  at  three  fixed 
rotational  angles  (9  = 0 , 45  , 90  ) tor  the  different 
attached  mass  ratios  Y = 0,  1,  5,  10,  to  show  the 
effect  of  mass  reinforcement  on  the  acoustic  pressure; 
Figu-e  3-4  was  made  at  two  rotational  angles 
(9  - 45  , 90=>  for  various  depths  of  the  attachment 
( q ^ = 0.  0.37,  1-85  , 3.7)  to  show  the  stiffness  efrect 
on  the  acoustic  pressure,  and  Figure  3.5  was  made  at 
three  fixed  rotational  angles  (0=0',  45J,  90  ) tor 
various  beam  sizes  to  show  both  the  influence  of  the  mass 
and  the  stiffness  on  the  radiated  acoustic  pressure 
(3)  Above  tiie  coincidence  frequency  at  (2  = 2 5:  Three 

figures,  i.e  , Figure  3.6,  3.7,  and  3.8  were  plotted 
in  the  same  manner  as  Figures  3.3,  3 4,  and  3.5, 
respectively . 

3.1.3  Power  Spectrum.  The  radiated  acoustic  power  7N) 
from  a point-excited  beam-reinforced  plate  was  computed  relative  to 
the  driving  power  (tO  of  an  infinite  isotropic  plate  excited  by  a 
point  force  at  the  origin,  versus  the  nondimensional  frequencies 
((2) . The  driving  power  (N  ) can  be  shown  to  be 
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and  therefore, 


Plate  (Stiffness  Only, 


BEAM 


Directivity  Function  of  Far-Field  Pressure  of  Beam-Reinforced 
Plate  (Beam,  £2  = 1.5) 


Directivity  Function  of  Far-Field  Pressure  of  Beam-Reinforced 
Plate  (Stiffness  Only,  ft  = 2.5) 
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Figure  39  presents  the  power  spectrum  oi  a plate  with  an 
attached  line  mass  for  five  mass  ratios  y * 0,  1 0,  2-7,  50,  10.0 
Figure  3.10  presents  the  power  spectrum  of  a plate  with  an  attached 
line  stiffness  for  five  various  depths  of  attachment  q = 0,  0 37, 
1.0,  1.85,  3.7.  Figure  3.11  presents  the  power  spectrum  of  a beam- 
reinforced  plate  for  various  beam  sizes. 


3 . 2 Numerical  Results  for  Beam-Reinforced  Plate  with  Two 

Additional  Spring-Mass  Systems 


3.2.1  Input  Data.  For  an  aluminum  infinite  plate  in 
contact  with  water,  numerical  results  are  presented  by  using  the  same 
physical  constants  given  in  Section  3-1  1, 

3.2  2 Far-Field  Pressure  Spectrum  and  Directivity  Function. 
The  far-field  pressure  in  the  present  problem  (see  Figure  3.3)  is 
shown  in  Equation  (2.74),  i e.. 
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NORMALIZED  FREQUENCY  = 


Normalized  Power  N/N  of  Beam-Reinforced  Plate  (Stiffness  Only) 
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Similar  to  the  expression  for  the  radiated  pressure  in 

Equation  (3.1),  the  bracketed  expression  in  Equation  (2.74)  is 

seen  to  be  the  radiated  far-field  pressure  from  a point-excited 

isotropic  plate.  The  nondimens ional  coefficient  C is  a pressure 

ms 

amplitude  correction  factor  for  a point-excited  beam-reinforcement 
plate  with  two  additional  spring-mass  systems  equally  located  from 
the  attached  beam.  This  coefficient  C is  much  more  complicated 
than  C in  Equation  (3.1)  due  to  the  couplings  of  some  additional 
relevant  factors  such  as  sprung  mass  (m^) , spring  stiffness  (k^)  and 
location  of  the  spring-mass  systems. 

If  the  excitation  frequency  of  the  point  force  approaches 
the  natural  frequency  of  the  spring-mass  system,  the  radiated 
pressure  approached  the  following  limit 
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where 


- [2I2°cos(ko<l  cos0  sin<J>)  - I^0-^0]  / [-I1°-I3°+{u1y-  u3Uysln40sm%} 

< I1°2+I10I30-2I2°2} 1 . (3.5) 

It  should  be  noted  that  as  u)  -*■  w^,  each  spring-mass  acts  like  a 
simple  support;  in  other  words,  the  transverse  displacement  under  the 
spring  is  zero,  but  rotation  of  the  plate  is  not  zero.  Thus,  the 
radiated  pressure  is  due  to  a point  excited  beam-reinforced  plate, 
simply  supported  at  x * +d. 

By  using  the  radiated  far-field  pressure  expression  in 
Equation  (2.74),  the  pressure  spectrum,  directivity  function,  and 
power  spectrum  were  computed  in  a similar  way  as  was  discussed  in 
Section  3.1.  It  should  be  noted  that  all  the  following  numerical 
computations  done  for  this  particular  problem  were  based  on  the 
fixed  dimension  of  the  attached  beam,  i.e.,  (y  “ 2.7,  - 1). 

(i)  The  far-field  pressure  spectra  at  zero  azimuthal  angle 
(<J>  • O’)  per  1 pound  force  at  1 yard  in  dB  relative  to  1 y bar-in 
were  computed  versus  normalized  frequencies  (ft)  for  several  different 
conditions: 
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Figures  3.12,  3.13,  and  3.14  show  the  pressure  spectra 
2 

of  fixed  spVung  mass  » 5.0)  and  spring  characteristic 

frequency  (aL/w^  ■ 0.1,  or  1.0, or  2.5)  with  the  variables  kQd  ■ tt/5, 


it,  3tt,  to  illustrate  the  influence  of  the  location  of  the  spring-mass 
systems  on  the  radiated  pressure.  Figures  3.15,  3.16  and  3.17 
show  the  pressure  spectra  of  a fixed  spring  mass  (m^/pQh  « 5.0) 
and  its  distance  (kQd  ■ tt/5 , or  n,  or  3tt)  with  variables  spring 
stiffness,  represented  by  * 0.1,  1.0,  2.5,  to  illustrate  the 

influence  of  the  spring  stiffness  k^  on  the  radiated  pressure. 

Figures  3.18,  3.19,  and  3.20  show  the  pressure  spectra  of  a fixed 

spring  characteristic  frequency  (co^/cu^  - 1.0)  and  the  distance  of  the 

2 

spring-mass  system  (kQd  - tt/5,  or  n,  or  3tt)  with  variables  m^/pQh  » 1.0, 
2.7,  5.0,  to  illustrate  the  effect  of  spring  mass  on  the  radiated 
pressure . 

fii)  The  directivity  function  for  the  radiated  far-fleld 
pressure  from  a point-excited  beam-reinforced  plate  with  two  additional 
spring-mass  systems  were  calculated  relative  to  the  far-field  pressure 
at  zero-azimuthal  angle  for  a point-excited  isotropic  plate,  versus 
the  azimuthal  angle  4>  at  three  rotational  angles  (0  - 0°,  45e, 

90°)  and  different  excitation  frequencies  (ft  - 0.01,  0.1,  0.5, 

1.0,  1.5,  2.0,  2.5): 
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Plots  of  the  directivity  function  were  made  for  Q = 0.5,  1.5  and  2.5 

only  and  are  presented  in  nineteen  figures. 

The  most  Interesting  directivity  plots  are  those  for  which 

the  excitation  frequency  (to)  coincides  with  the  characteristic 

frequency  (uj^)  of  the  spring-mass  system.  Three  figures  of  this 

case,  i.e..  Figures  3.21,  3-22,  and  3.23  are  plotted  at  frequency 

(ui./u  = w/u ) =0.5,  or  1-5  or  2.5)  to  show  the  location  effect 

1 o o 

of  the  spring-mass  system  on  the  resulting  pressure. 

Five  figures  of  the  directivity  function  were  plotted  for  a 

frequency  below  the  coincidence  frequency,  i.e.,  at  = 0 5,  which 

are  shown  in  Figures  3.24  to  3.28.  Figure  3.24  shows  the  directivity 

function  at  fixed  sprung  mass  and  characteristic  frequency 
2 

(m^/p^h  =2  7,  = ^ f°r  different  locations  of  the  spring- 

mass  system  (kod  * 0,  rr/ 5 , n) . Figure  3-25  and  3.26  show  the 

2 

directivity  function  at  fixed  spring  mass  (m^/p^h  *>  2.7)  and  its 

location  (k^d  = tt/5  or  3^)  for  different  spring  characteristic 

frequencies  * 0,  0.1,  1.0,  2.5).  Figure  3.27  and  3.28 

show  the  directivity  function  at  the  fixed  spring-maws  location 

(k  d = W5  or  3tt)  and  its  characteristic  frequency  (oo,  /co  = 1.0) 
o 1 o 

2 

at  different  sprung  masses  (m./Poh  « 0,  2-7,  5.0). 

Six  figures  of  directivity  functions  were  plotted  for  a 

frequency  above  the  coincidence  frequency,  i.e.,  at  ft  * 1.5. 

Figures  3.29  and  3.30  show  the  directivity  function  at  fixed  sprung 
2 

mass  (m^/pQh  = 2.7)  and  its  characteristic  frequency  (co^/ooq  - 1.0, 
or  2.5)  for  different  locations  (kQd  ■ 0,  tt/5,  tt , or  3n).  Figures 


Fairing  Model  (ft 


Figure  3.22  Directivity  Function  of  Far-Field  Pressure  of  an  Acoustic 
Fairing  Model  (ft  = 1.5,  W./u  = 1.5,  m /p  h2  = any) 


Figure  3.24  Directivity  Function  of  Far-Field  Pressure  of  an  Acoustic 
Fairing  Model  (Q  * 0.5,  u>. /to  = 1.0,  m /p  = 2.7) 


Figure  3.25  Directivity  Function  of  Far-Field  Pressure  of  an  Acoustic 
Fairing  Model  (ft  = 0.5,  k d = tt/5,  m../p  h2  = 2.7) 


106 


Figure  3.26  Directivity  Function  of  Far-Field  Pressure  of  an  Acoustic 
Fairing  Model  (fi  = 0.5,  k d = 3ir,  m.  /p  = 2.7) 
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Figure  3.29  Directivity  Function  of  Far-Field  Pressure  of  an  Acoustic 
Fairing  Model  (Q  = 1.5,  oo./u)  = 1.0,  m /p  h2  = 2.7) 


ill 


3.31  and  3.33  show  the  directivity  function  at  fixed  sprung  mass 
2 

(m./p  h = 2.7)  and  its  location  (k  d = tt/5,  or  3tt)  for  different 
i o o 

spring  characteristic  frequencies  = 0,  01,  1.0,  2.5). 

Figures  3.32  and  3 34  show  the  directivity  function  at  the  fixed 

spring  characteristic  frequency  = 10)  and  its  location 

2 

(k^d  = rr / 5 , or  3TT)for  different  sprung  masses  (m^/p^h  = 0,  2.7, 

5.0). 

Five  figures  of  directivity  function  were  similarly  plotted 
for  a frequency  above  the  coincidence  frequency,  i e , at  57  * 2.5, 
similar  to  those  plotted  for  51  “ 1.5.  Figure  3.35  was  plotted 
for  a fixed  sprung  mass  and  Its  characteristic  frequency  for  different 
locations;  Figures  3.36  and  3.38  were  plotted  for  a fixed  sprung 
mass  and  its  location  for  different  spring  characteristic  frequencies; 
Figure  3.37  and  3.39  were  plotted  for  a fixed  spring  characteristic 
frequency  and  its  location  for  different  sprung  masses. 

3.2.3  Power  Spectrum.  The  radiated  acoustic  power  (N)  frcm 
a point-excited  beam-reinforced  plate  with  two  additional  spring- 
mass  systems  was  computed  relative  to  the  driving  power  (N^)  of 
an  infinite  isotropic  plate  excited  by  a point  force  at  the  origin, 
versus  the  nondimensional  frequencies  (57)  at  different  conditions. 


Function  of  Far-Field  Pressure  of  an  Acoustic 


Directivity  Function  of  Far-Field  Pressure  of  an  Acoustic 
Fairing  Model  (SI  = 1.5,  k d = 3ir,  m1 /p  = 2.7) 


Figure  3.34  Directivity  Function  of  Far-Field  Pressure  of  an  Acoustic 
Fairing  Model  (51  = 1.5,  k d = 3n,  w /oj  = 1.0) 
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Function  of  Far-Field  Pressure  of  an  Acous 


Function  of  Far-Field  Pressure  of  an  Acous 


5 2 = U mm  n = °<»/ 


Figure  3.38  Directivity  Function  of  Far-Field  Pressure  of  an  Acoustic 
Fairing  Model  (ft  = 2.5,  k d = 3^,  m./p  = 2 
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Eight  figures  were  plotted  for  the  power  spectrum.  Figure 

3.40  was  made  at  fixed  spring  characteristic  frequency  to^/u^  * 0*1 

2 

and  sprung  mass  m, /p  h =2.7  for  different  locations  k d = tt/5,  tt,  3it. 
r ° 1 o o 

Figure  3.40  is  also  valid  for  the  case  of  fixed  spring  characteristic 

frequency  ■ 1 and  location  kQd  = tt  for  different  sprung 

2 

mass  ratios  m, /p  h =1.0,  2.7,  5.0.  Figure  3.41  and  3.42  were  made 
1 o 

at  fixed  spring  characteristic  frequency  = 1.0  or  2.5)  and 

2 

sprung  mass  (m^/p  h = 2.7)  for  different  locations  (kQd  = tt/5,  3tt)  ; 

2 

Figure  3.43,  3.44,  and  3.45  were  made  at  fixed  sprung  mass  (m^/PQh  = 

2.7)  and  its  location  (k  d = tt/5  or  tt,  or  3tt)  for  different  spring 

o 

characteristic  frequencies  (w^/u)q  = 0.1,  1.0,  2.5);  Figure  3.46  and 

3.47  were  made  at  fixed  spring  characteristic  frequency  (co^/u^  = 1.0 

2 

or  2.5)  for  different  sprung  masses  (m^/p^h  = 1.0,  2.7,  5.0). 


NORMALIZED  FREQUENCY  SI  = w /u 


Normalized  Power  N/N  of  an  Acoustic  Fairing 
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CHAPTER  IV 


DISCUSSIONS  OF  THE  RESULTS  AND  CONCLUSIONS 
4 . 1 Introduction 

The  initial  task  in  the  analysis  of  a vibrating  marine 
structure  is  to  initially  make  some  simplifying  assumptions  about 
its  motion  and  then  idealize  the  structure  so  that  the  analysis  of  its 
motion  can  be  made  feasible.  In  this  thesis,  an  idealized  beam- 
reinforced  plate  will  be  conceived  as  a ship's  hull  reinforced  by 
a framing  element  (beam) , and  a time-harmonic  point  force  represents 
the  action  of  a vibrating  machine. 

Vibrating  machinery  is  usually  supported  by  frames  of  a 
ship's  hull  due  to  structural  design  considerations.  Since  the 
reduction  in  the  radiated  noise,  when  a vibrating  machine  is 
attached  to  a reinforcing  beam  of  a hull-plating  has  not  been 
analytically  proven,  the  first  major  task  in  this  thesis  was  to  study 
the  acoustic  fields  produced  by  a point  force  excitation  of  an 
infinite  beam  which  is  welded  to  an  infinite  plate.  The  attached 
beam  is  considered  to  act  as  a structural  discontinuity  of  one  of  the 
following  three  types:  (1)  a mass  discontinuity  which  will  normally 

act  as  a line  mass  Impedance,  (2)  an  elastic  discontinuity  which  acts 
as  a path  of  energy  propagation,  and  (3)  a beam-like  discontinuity 
which  acts  as  an  elastic  as  well  as  mass  discontinuity.  This  problem 
constitutes  an  initial  step  leading  to  a better  understanding  of  the 
role  a structural  discontinuity  plays  in  an  infinite  reinforced 
plate.  Besides,  a mathematical  procedure  for  solving  the  sound 
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radiation  of  an  infinite,  reinforced  plate  developed  for  this 
problem  can  be  utilized  to  obtain  the  sound  radiation  from  other 
plates  with  more  complex  structural  discontinuities.  Numerical 
results  of  this  problem  were  plotted  and  shown  from  Figures  3.1  to 
Figure  3.11.  These  will  be  subsequently  discussed  in  detail  in 
Section  4.2.  A second  major  task  of  this  thesis  was  to  devise  a 
simple  mathematical  model  of  acoustic  fairing,  i.e.,  by  attaching  two 
line  spring-mass  systems  to  the  plate  as  shown  in  Figure  2.2,  and 
to  study  their  interaction  with  the  existing  attached  beam.  The 
sound  radiation  from  an  acoustically  faired  plate  was  investigated 
by  varying  the  various  parameters  that  influence  it,  such  as  the 
mass  and  the  spring  stiffness  as  well  as  the  location  of  the  systems 
with  respect  to  the  reinforcing  beam.  This  mathematical  model, 
though  simple,  yields  a better  understanding  of  the  influence  of  the 
physical  parameters  on  the  application  of  acoustic  fairing.  The 
numerical  results  obtained  for  this  second  problem  were  plotted 
from  Figures  3.12  to  3.47.  The  detailed  discussions  of  this  problem 
are  given  in  Section  4.3. 

4.2  Far-Field  Acoustic  Radiation  from  Point  Excited  Beam-Reinforced 

Plate 

The  influence  of  structural  discontinuities  on  the  radiated 
sound  from  plates  attached  to  such  discontinuities  was  assessed 
mathematically  in  Section  2.3  and  many  numerical  examples  were 
presented  in  Section  3.1. 

The  influence  of  a mass,  a stiffness,  or  a beam  on  the 
radiated  sound  can  be  best  discussed  in  the  two  frequency  ranges: 
below  and  above  the  coincidence  frequency. 


n 


i 
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4.2.1  Excitation  Frequency  Below  the  Coincidence  Frequency 

(ft  < 1) . Since  the  frequency  of  excitation  is  below  the  coincidence 

frequency,  one  does  not  expect  a peak  in  the  radiated  far-field 

pressure.  Thus,  maximum  pressure  is  observed  at  points  directly 

across  the  point  of  excitation,  i.e.,  <p  = 0°. 

(i)  Influence  of  Mass:  The  influence  of  an  infinitely  long 

uniformly  distributed  mass  on  the  radiated  acoustic  pressure  is 

exhibited  by  the  factor  C of  Equation  (3.1).  If  one  sets  U =0 

ms  y 

in  this  expression  in  Equation  (3.1)  and  obtains  an  approximate  value 

for  I^°  for  low  frequencies,  it  was  shown  that  this  product 

0 8 

increases  with  the  frequency  as  (yft  ' ) . This  indicates  that  the 

factor  will  decrease  in  value  if  the  frequency  ratio  0 or  the 

mass  ratio  y increases.  The  reduction  of  the  radiated  pressure  due 

to  the  addition  of  an  infinite  line  mass  is  clearly  shown  in  Figure  3.1 

A mass  ratio  of  y = 10,  which  is  equivalent  to  a mass  whose  cross- 

2 

sectional  area  is  equal  to  3.7h  , reduces  the  radiated  pressure  by 
5 dB  at  best,  for  frequencies  ft  < 1. 

The  influence  of  the  distributed  mass  on  the  directivity 
function  at  frequencies  below  the  coincidence  frequency  can  be 
expected  to  be  minimal.  This  is  shown  clearly  by  the  directivity 
plots  at  0 = 0°  and  90’  in  Figure  3.2. 

The  influence  of  the  mass  on  the  radiated  power  is  shown  in 
Figure  3.9.  It  is  apparent  that  one  can  attain  a reduction  in  the 
radiated  power  equivalent  to  that  of  the  radiated  pressure,  because  of 
the  omnidirectionality  of  the  directivity  function.  For  a mass  ratio 
Y - 10,  the  radiated  pcver  was  reduced  by  as  much  as  4.5  dB  for 
frequencies  ft  < 1.  The  reason  for  reduction  in  the  radiated  power  can 
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be  attributed  to  the  fact  that  the  mass  requires  additional  power  to 
accelerate;  hence,  a smaller  portion  of  the  input  power  is  pumped  to 
the  acoustic  medium. 

(ii)  Influence  of  Stiffness:  The  influence  ci  an  infinitely 

long  line  stiffness  on  the  radiated  pressure  is  exhibited  by  the 
factor  Cmg  in  Equation  (3-1).  If  one  sets  y =>  Cr  in  Equation  (3  1)  and 
obtains  an  approximate  expression  for  the  products  I^°  for  low 

frequencies,  it  was  shown  that  the  product  increases  with  the  frequency 

2 8 — 
as  (U  12  ' ) , This  indicates  that  the  factor  C will  decrease  as 
y ms 

follows: 


mS  1 + (...)i  U (2^  ® sin^0  sin\ 

y 


At  the  angle  of  observation  0=0°,  the  factor  C is  not  influenced 

ms 

by  the  stiffness  parameter  and  is  equal  to  unity.  Thus,  che 
pressure  spectrum  in  Figure  31  does  not  exhibit  the  influence  of  the 
stiffness . 

The  influence  of  the  line  stiffness  on  the  directivity 
function  is  minimal  for  frequencies  below  the  coincidence  frequency, 
as  shown  by  the  directivity  plots  at  0 = 0°  and  90°  in  Figure  3.2. 

The  influence  of  the  line  stiffness  on  the  radiated  power 
is  shown  in  Figure  3.10.  It  is  apparent  that  a reduction  of  radiated 
power  by  2.5  dB  was  attained  when  a line  stiffness  having  a depth  ratio 
of  q^  * 3.7  (or  a moment  of  Inertia  ratio  of  45)  is  attached  to  the 
plate.  The  reduction  in  the  acoustic  radiated  power  can  be  attributed 
to  the  additional  power  required  by  the  elastic  (mechanical)  power 
of  the  infinite  elastic  stiffener  and  hence  less  of  the  input  power 
is  available  for  the  acoustic  medium. 
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(iil)  Influence  of  a Beam:  The  influence  of  an  elastic 

beam  welded  to  an  infinite  plate  on  the  radiated  pressure  is  exhibited 
in  Equation  (3.1).  The  complicated  interaction  of  the  mass  of  the  beam 
and  its  stiffness  is  shown  primarily  by  the  factor  C of  Equation 

mS 

(3.1).  However,  below  the  coincidence  frequency  (I)  < 1),  the 

0 8 

influence  of  the  mass  is  proportional  to  ft  ' while  that  of  the 

2 g 

stiffness  is  proportional  to  ft  ‘ , which  indicates  that  the  stiffness 
of  a reinforcing  beam  of  reasonable  dimensions  is  overshadowed  by 
the  influence  of  its  mass. 

Since  the  stiffness  of  a beam  has  no  influence  on  the  radiated 

sound  spectrum  at  4»  = 0°,  the  only  factor  that  dominates  the  sound 

radiation  at  this  observation  point  is  due  to  its  mass,  as  is  clearly 

shown  in  Figure  3.1-  Thus,  the  mass  of  the  reinforcing  beam  does 

reduce  the  radiated  pressure  by  5 dB  for  a mass  ratio  y = 10. 

The  influence  of  the  beam  on  the  directivity  function  can  be 

expected  to  be  minimal  below  the  coincidence  frequency  as  exhibited 

by  the  directivity  plots  at  6 » 0°  and  0 ■ 90°  in  Figure  3.2. 

The  influence  of  the  reinforcing  beam  on  the  radiated  power 

is  shown  in  Figure  3.11.  The  reduction  of  radiated  sound  power  can 

reach  as  much  as  5 dB  for  a beam  of  reasonable  dimensions.  Below 

the  coincidence  frequency,  the  radiated  power  reduction  attained  by 

the  attachment  of  a beam  seems  to  be  dominated  by  its  mass  rather  than 

by  its  stiffness,  as  indicated  by  the  dependence  of  the  former  on  the 

0.8  28 

frequency  as  ft  ' and  the  latter  as  ft  ’ for  ft  < 1.  It  is  interesting 
to  note  that  Howe  and  Heckl  [Reference  13,  Equation  (38)]  make  the 
same  point  when  they  discuss  the  influence  of  random  fluctuations  of 
plate  mass  and  on  the  radiated  sound  from  the  plate  when  excited  by  a 
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point  force.  Thus,  the  major  part  of  the  5 dB  noise  power  reduction 
due  to  the  attachment  of  a beam  can  be  attributed  to  its  mass 
(4.5  dB)  than  its  stiffness  (2.5  dB)  . 

An  estimate  of  the  radiated  power  from  an  excited  beam- 
reinforced  plate  can  be  obtained  by  assuming  that  the  beam  is  welded 
to  the  plate  at  only  one  point,  the  point  of  excitation.  This 
obviously  neglects  the  complex  interaction  of  the  impedance  of  the 
beam  along  its  length  with  that  of  the  plate  along  an  Infinite  line. 

It  should  be  expected  that  the  coupling  of  the  beam  to  the  plate  is 
minimal  for  shallow  and  strong  for  deep  beams.  Denoting  the 
characteristic  impedance  of  an  unloaded  plate  and  beam  by  Zp  and  , 
then  the  portion  of  the  force  available  for  the  excitation  of  the  plate 
is  given  by  Z / (Z  + Z,  ).  Using  this  force  as  an  excitation  force  on 

p p D 

an  acoustically  loaded  isotropic  plate,  the  ratio  of  the  power 
radiated  from  the  beam-plate  to  that  radiated  by  an  isotropic  plate 
below  the  coincidence  frequency  can  be  expressed  by  the  square  of 
this  ratio,  see  Skudryzk  [15,  Eq.  (12.209)].  Thus, 

^plate-beam 

N , - 
plate 

Expressions  for  the  driving  point  impedance  of  an  infinite  plate  and 
an  infinite  beam  are  given  in  Skudryzk  [15,  Equations  (6.90)  and 
(9.45))  as  follows: 

Zp  - 8i/mD  , (4.2) 

Zb  - 2 (1+i)  M/U 


Z + Z, 


(4  1) 


(4.3) 
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The  estimated  power  ratio  as  expressed  in  Equation  (A.l)  is  expressed 
by: 

Power  Ratio  = /3  — q (^)  v^T  . (A.  A) 

z c z h 

P 

The  expression  for  the  power  ratio  in  Equation  (A. A)  is  proportional 
to  the  beam-plate  depth  ratio  (q ^ ) and  the  beam  breadth-plate  depth 
ratio  (b/h)  as  well  as  the  frequency  ratio  Si(“  ) . The  estimate 

of  power  ratio  given  in  Equation  (A. A)  can  be  compared  to  the  exact 
ratio  taken  from  Figure  3.11  as  follows  for  the  case  of  b = h. 

Table  A.l 


Power  Ratios  of  Beam-Reinforced  Plate 
(Below  Coincidence  Frequency) 


Depth  Ratio 
q2  » hb/h 

Stiffness 

Ratio 

I.  /bl 
b 

Power  Ratio  dB 

Exact , 

SI  - 0.5 

Estimate , 
Si  = 0.5 

Exact, 

Si  = 1.0 

Estimate , 
Si  = 1.0 

1 

1 

0.6 

l.A 

1.0 

3.0 

1.8 

'5.6 

1.0 

3.3 

2.2 

A.  A 

3.7 

A5 

2.0 

6.6 

5.0 

8.6 

The  estimate  of  radiated  power  for  beams  which  are  shallow  is  within 
2 dB  of  the  exact  solution.  The  plate  impedance  did  not  account  for 
the  fluid  mass  loading  on  the  plate. 

A. 2. 2 Excitation  Frequency  Above  the  Coincidence  Frequency 
(Si  > 1) . The  radiated  acoustic  pressure  above  the  coincidence 


frequency  has  a peak  at  a specific  azimuthal  angle.  The  peak  in  the 
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radiated  pressure  from  a point  excited  infinite  isotropic  plate 
occurs  at  the  angle  0^: 


sin  <f> 

c 


— for  ft  > 1 


(4.5) 


The  maximum  pressure  at  the  coincidence  angle  0^  for  an  isotropic 
plate  can  be  obtained  by  substituting  Equation  (4.5)  into  Equation 
(3.1)  and  setting  Cmg  = 1,  resulting  in  the  following  expression: 


max 


-i  k F 
o o 

2ttR 


. 1/2  ik  R 
,,  lx  O 

(1  " ft}  6 


(4.6) 


When  the  plate  is  reinforced  by  an  attached  mass,  stiffness  or  beam, 
the  peak  in  the  radiated  pressure  shifts  in  azimuthal  angle  and  the 
maximum  pressure  changes  in  magnitude.  There  seems  to  be  no  closed 
form  formulas  for  the  prediction  of  the  peaks  or  the  azimuthal 
coincidence  angles  for  reinforced  plates. 

(i)  Influence  of  Mass:  The  influence  of  an  infinite  line 

mass  on  the  radiated  pressure  is  exhibited  by  the  factor  C^s  in 
Equation  (3.1).  If  one  sets  * 0 and  obtains  an  approximate 
expression  for  the  factor  1^°  y^y,  an  approximate  expression  for  C^g 
valid  for  ft  > 1 and  — 0°  , can  be  obtained  as  follows: 


C 

ms 


1 - u1  (1  + 0.09yft1/2)  - i(0.09yft1/2) 


(4.7) 


Examination  of  the  factor  C , given  in  Equation  (4.7),  indicates  tht 
it  decreases  monotonically  with  mass  ratio  and  excitation  frequency. 
This  is  clearly  shown  by  the  numerical  results  presented  for  the 
radiated  pressure  spectrum  at  0 * 0®,  shown  in  Figure  3.1.  To  obtain 
a crude  estimate  of  the  reduction  in  the  radiated  pressure  at  0 « 0°, 
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one  can  resort  to  physical  arguments.  If  the  attached  mass  is  spread 
over  a wavelength  X and  a new  plate-mass  surface  density  is  computed 
and  used  in  the  expression  for  the  radiation  from  an  infinite 
isotropic  plate,  one  may  get  a crude  estimate  for  the  pressure.  Thus, 


the  new  surface  mass  density  in  becomes: 

2 


M 


Pq2h' 


m = m + y=Ph+  y 


q2h 

m ( 1 + -~)  , 


where  X is  the  bending  wavelength  in  the  plate  given  by 


, 2tt  1 4 

X = — a,  a 

/jT~ 

o 


h2E 


12p (1-V  ) 


2 lm~ 

°o  = c Y D 


/12  c2 
h c 
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(4.8) 


Thus,  the  new  mass  density  m becomes: 


m = m(l  + — q„  ) 

7T  c 2 

P 


(4.9) 


The  radiated  pressure  from  an  isotropic  plate  p^  can  be  written  as 
follows : 


-ik  F 
o o 


ik  R 

o 

e cos 


o 2ttR  ik  mcostj) 

1 (1  - f2  sin  $) 


(4.10) 


Thus,  the  ratio  of  the  radiated  pressure  of  a plate  with  an  attached 
mass  p^  to  that  of  an  isotropic  plate,  p^,  computed  at  $ » O'  can 
be  written  as: 


2~2,_ 1/2 


(4.11) 


1 '39 


The  estimate  given  in  Equation  (4.11)  can  be  compared  to  the  exact 
solution  shown  in  Figure  3.1  in  the  following  table. 

Table  4.2 

Pressure  Ratios  of  Mass-Reinforced  Plate 


p / p in  dB 

o o 

Mass  Ratio 
Y 

Exact, 
ft  = 1.5 

Estimate,  Exact, 

ft  - 1.5  ft  - 2.5 

Estimate , 
ft  - 2.5 

1 

0.8 

0.6 

1.0  0.8 

10 

6.5 

4.6 

8.2  5.8 

Table  4.2 

shows  that  the 

reduction  of  the 

radiated  pressure  spectrum 

at  $ • 0° 

can  reach  8 dB 

for  a mass  ratio 

of  y » 10.  The  table  also 

shows  that  the  crude  estimate  given  by  Equation  (4.11)  does  not 
exceed  2.5  dB  at  the  upper  limit  of  mass  and  frequency  ratios. 

The  directivity  function  for  the  radiated  pressure  has  a peak 
at  the  same  azimuthal  angle.  Since  there  is  no  way  one  can  predict 
this  angle  analytically,  a crude  estimate  of  the  coincidence  angle 
can  be  made.  It  is  clear  that  the  mass  would  not  alter  the  coincidence 
angle  at  the  plane  9 » 0J;  hence,  the  coincidence  angle  ^ at  9 • 0° 
plane  is  given  by  Equation  (4.5).  The  directivity  function  will  be 
most  affected  near  or  at  the  plane  8 » 90'.  This  is  exhibited  by 
Figures  3.3  and  3.6.  If  the  attached  mass  is  spread  over  a wavelength 
and  the  resulting  mass  density  added  to  that  of  the  plate,  one  can 
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obtain  a new  coincidence  frequency  co^  as  follows: 


(4.12) 


Using  the  new  coincidence  frequency  in  the  expression  for  the 
critical  angle  <p^  in  Equation  (4.5)  results  in 

sin  4.  =—  [1  + 0 q_  ]1/4  . (4.13) 

C Al  n cp  2 

Equation  (4.13)  can  be  used  to  estimate  the  new  coincidence 
angle  , which  is  then  compared  to  the  exact  solution  as  exhibited 
in  Figures  3.4  and  3.6.  The  comparison  is  shown  in  Table  4.3. 


Table  4.3 


New  Coincidence  Angles  for 

Mas s-Re info reed 

Plate 

Mass  Ratio 

*c 

for  ft  = 1.5 

*c 

for  ft 

= 2.5 

Y 

Exact 

Estimate 

Exact 

Estimate 

0 

55° 

55° 

39° 

39“ 

1 

55° 

56° 

40° 

40“ 

5 

60° 

62° 

45“ 

44° 

10 

66° 

70“ 

52° 

•ts 

00 

o 

Comparison  of  the  exact  and  estimated  coincidence  angles  show  a good 

agreement  even  at  the  high  frequenoy  and  mass  ratios.  In  general,  an 

attached  mass  increases  the  magnitude  of  the  coincidence  angle  4>  . 

c 


'"W; 


v 
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The  magnitude  of  the  pressure  peak  at  the  new  coincidence 

angles  decreases  as  y increases  at  the  plane  9 = 0C  . This  is  due  to 

the  fact  that  at  0 = 0°,  the  expression  for  the  ratio  of  the 

directivity  functions  of  a plate  with  an  attached  mass  to  that  of  an 

isotropic  plate  is  the  same  as  that  for  <p  = 0°  given  in  Equation  (4.7). 

.n  fact,  comparison  of  these  ratios  at  peak  pressures  gave  the  same 

dB  levels  as  presented  in  Table  4.2.  The  magnitude  of  the  pressure 

peak  at  9 = 90°  cannot  be  predicted  as  easily  as  that  at  9 =0°. 

It  can  be  seen  from  Figures  3.3  and  3.6  that  the  pressure  peak  for 

a plate  with  an  attached  mass  can  be  higher  or  lower  than  that  of  an 

isotropic  plate.  It  is  apparent  that  the  coefficient  C becomes 

ms 

very  large  when  <p  approaches 

The  radiated  power  from  a plate  with  an  attached  mass  is 

shown  in  Figure  3.9.  The  radiated  power  decreases  with  increasing 

attached-mass  ratio  y,  the  reduction  of  the  radiated  power  can  be  as 

high  as  5 dB  at  fi  = 1.5  and  3 dB  at  0 - 2.5  and  y = 10. 

(ii)  Influence  of  Stiffness:  The  influence  of  an  infinite 

line  stiffness  on  the  radiated  pressure  is  exhibited  by  the  factor 

C in  Equation  (3.1).  If  one  sets  y = 0,  and  <)>  = 0“ , then  C =1.0, 
ms  ms 

which  indicates  that  the  pressure  spectrum  at  <f>  = 0°  for  a plate 
with  an  attached  line  stiffness,  is  the  same  as  that  for  an  isotropic 
plate . 

The  directivity  function  for  the  radiated  pressure  of  a plate 
with  an  attached  infinite  line  stiffness  has  a peak  at  some  azimuthal 
angle.  Since  there  is  no  way  to  predict  this  angle  analytically,  a 
crude  estimate  of  the  coincidence  angle  can  be  made.  It  is  obvious 
that  the  line  stiffness  does  not  alter  the  coincidence  angle  at  the 
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plane  6 « 0 , because  the  directivity  function  at  0 = 0°  is  the  same 
as  that  for  an  isotropic  plate.  The  directivity  function  will  be 
changed  considerably  near  the  plane  0 = 90°.  This  effect  is  exhibited 
in  Figures  3.4  and  3.7.  If  the  attached  stiffness  is  spread  over  a 
wavelength  X and  a new  plate  with  a new  stiffness  is  computed,  the 
ratio  of  the  stiffness  of  the  reinforced  plate  D to  that  of  an 
isotropic  plate  D becomes: 


D 

D 


1 + q 22C 


3£(l+q22) 

+ I H 


(4.14) 


where 

5 * q2h/x  m 7~  ^ • 

p 

Using  the  ratio  of  the  new  stiffness  D to  the  isotropic  stiffness  D, 
a new  coincidence  frequency  can  be  computed  as  follows: 


Vf-.VT 


(4.15) 


Using  the  new  expression  for  tuo>  a crude  estimate  of  the  coincidence 


angle  <j>  is  obtained  as  follows: 

, 35(1  + q 2) 

sin  <P  " — [1  + q £ + - 

c /ft  2 


-1/4 


1H 


(4.16) 


Equation  (4.16)  was  used  to  compute  an  estimate  for  the  coincidence 
angles  which  are  then  compared  with  the  exact  values  as  presented 

in  Figures  3.4  and  3.7.  The  comparison  of  the  exact  to  the  estimated 
value  is  presented  in  Table  4.4. 
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Table  4.4 


New  Coincidence  Angles  for  Stiffness-Reinforced  Plate 


Stiffness  Ratio 

<P  for  ft  = 1.5 
c 

<(>  for  ft  * 2.5 
c 

q2 

Exact  Estimate 

Exact  Estimate 

0 

55° 

55° 

39  o 

39° 

0.37 

50° 

50° 

38° 

36° 

1.85 

25° 

32° 

19° 

23° 

3.7 

17° 

20° 

12° 

15s 

Table  4.4  shows  a good  agreement  between  the  estimated  coincidence 
angles  and  those  obtained  by  exact  methods.  It  can  be  seen  that  an 
attached  stiffness  decreases  the  magnitude  of  the  coincidence  angle  <}>^. 

The  peak  pressure  in  the  directivity  function  generally 
decreases  with  increasing  stiffness.  This  can  be  seen  from  Figures  3.4 
and  3.7.  The  reduction  in  the  peak  pressure  reaches  6 dB  at  ft  = 1.5 
and  11  dB  at  ft  * 2.5  for  ” 3.7. 

The  radiated  power  from  a plate  with  an  attached  infinite 
line  stiffness  was  computed  and  shown  in  Figure  3.10.  It  can  be 
seen  that  the  radiated  power  decreases  with  increasing  stiffness  of 
the  attached  line  stiffness.  The  reduction  in  radiated  power 
reaches  4 dB  at  ft  - 1.5  and  2.5  and  " 3.7. 

(iii)  Influence  of  a Beam:  The  influence  of  an  attached 

infinite  beam  on  the  radiated  pressure  is  exhibited  by  the  factor  C 
of  Equation  (3.1).  Since  the  pressure  spectrum  at  4>  = 0°  is  only 


■«- 
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influenced  by  the  mass  of  the  beam,  nothing  new  can  be  added  to  the 
discussion  already  presented  under  Item  (i)  above,  as  can  be  seen  in 
Figure  3.1. 


The  directivity  function  at  the  plane  9 = 90°  of  the  pressure 
has  a peak  at  some  azimuthal  angle  which  cannot  be  predicted  by 
analytical  means.  To  obtain  a crude  estimate  for  the  new  coincidence 
angle,  both  the  increased  stiffness  and  mass  can  be  distributed 
over  a wavelength  A.  The  new  coincidence  frequency  becomes 


(4.17) 


Thus,  the  new  coincidence  angle  0^  can  be  estimated  by  the  equation 


sin  0 = — 

c — 


i_  r i+i 

v®\  2 31 

L1  + q^  + — 


U/4 


311+q2  ) 


(4.18) 


where  £ was  defined  in  Equation  (4.14). 

The  expression  in  Equation  (4.18)  predicts  an  estimate  of  the 

new  coincidence  angles  0 , which  is  tabulated  in  Table  4.5.  These 

estimates  are  compared  to  the  exact  values  of  0 obtained  from 

c 

Figures  3.5  and  3.8.  A good  agreement  between  estimated  and  exact 

values  of  0 is  shown  in  Table  4.5.  It  can  be  seen  from  the 
c 

results  presented  in  Table  4.5  as  well  as  the  expression  for  0^  in 
Equation  (4.18)  that  the  stiffness  is  more  dominant  than  the  mass  of 
the  beam.  The  dominance  of  the  stiffness  over  the  mass  may  be 
better  understood  in  the  light  of  Equations  (4.17)  and  (4.18), 
since  the  mass  increases  linearly  with  the  depth  ratio  q2  but  the 
stiffness  increases  with  the  third  power  of  q2 . 
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Table  4.5 


New  Coincidence  Angles  for  Beam-Reinforced  Plate 


*c 

for  ft  ■=  1.5 

< p for  ft 
c 

= 2.5 

Beam 

Exact 

Estimate 

Exact 

Estimate 

Y * 0, 

q2  - 0 

55° 

55° 

39° 

39° 

Y - 2.7, 

q2  * 1.0 

40° 

45° 

30° 

32° 

Y * 5 , 

q2  = 1.85 

32° 

35° 

23° 

25° 

Y * 10, 

q2  = 3.7 

25° 

24® 

18° 

18° 

The  peak  pressure  at 

the  coincidence 

angle  $ cannot 

be 

predicted 

easily. 

One 

can  observe  that  the  peak  pressure,  in 

general , 

decreases 

with  increasing  mass 

and  stiffness. 

The  reduction 

in  peak 

pressure 

at  <(>  for 
c 

the 

plane  0 

= 0°  reaches  6. 

5 dB  at  ft  = 1.5 

9 

8 dB  at  Jl 

- 2.5  for 

q2 

= 3.7. 

Similar  reductions  in  the  peak 

pressure 

at  for 

c 

the  plane  9 

■ 90°  reaches  6 

dB  at  ft  « 1.5 

and  12  dB 

at  fi  - 2. 

5 for  q2  « 

3. 

7. 

The  radiated  power  from  a plate  with  an  attached  beam  was 
computed  and  shown  in  Figure  3.11.  In  general,  the  radiated  power 
decreases  with  increasing  stiffness  and  mass  of  a beam.  Using  the 
expression  for  the  estimated  radiated  power  given  in  Equation  (4.4), 
the  following  results  were  obtained: 
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Table  4.6 


Power  Ratio  of  B :am-Reinf orced  Plate 
(Above  Coincidence  Frequency) 


Beam 

Power  Ratio 

(dB) 

ft  - 

1.5 

ft  - 

2.5 

Exact 

Estimate 

Exact 

Estimate 

Y - 2.7,  q2  * 1.0 

2.0 

2.3 

2.7 

5.0 

Y * 5.0  q2  ^1.85 

3.7 

5.3 

4.6 

6.5 

Y - 10.0,  q2  - 3.7 

6 .6 

9.7 

7.6 

13.6 

The  agreement  between  estimated  and  exact  computation  of  the  radiated 
power  is  within  3 dH  for  all  beams  at  ft  = 1.5  and  within  6 dB  for  all 
beams  at  ft  - 2.5.  The  deeper  the  beam,  the  worse  the  estimate 
became.  However,  the  reduction  in  radiated  power  was  shown  to  exceed 
6 dB  at  « 1.5  and  8 dB  at  ft  « 2.5  for  the  deep  beam  = 3.7. 

4*2.3  Summary . In  this  section,  a summary  of  the  discussion 
in  Section  4.2  is  presented. 

(i)  Influence  of  an  Attached  Mass:  In  general,  the 

attached  mass  reduces  the  radiated  pressure  spe.trum  at  $ * 0° 
above  and  below  the  coincidence  frequency  by  as  much  as  8 dB  for  a 
moderately  heavy  mass  for  (1  £ 2.5. 

The  directivity  function  is  slightly  influenced  by  an 
attached  mass  below  the  coincidence  frequency.  The  directivity 
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function  is  highly  influenced  by  the  attached  mass  only  near  9 = 90°, 
where  the  coincidence  angle  increased  by  as  much  as  13°  at  the  upper 
limit  of  SI  » 2.5.  The  peak  of  pressure  could  increase  up  to  2 dB. 

The  radiated  power  decreases  with  increasing  mass  of  the 
attached  line  mass,  reaching  as  high  as  5 dB  for  a moderately  heavy 
mass  for  SI  <_  2 . 5 

(ii)  Influence  of  an  Attached  Stiffness:  The  attached 

stiffness  has  no  influence  on  the  radiated  pressure  spectrum  at 
4>  = 0s  above  and  below  the  coincidence  frequency. 

The  directivity  function  is  slightly  influenced  by  an 
attached  stiffness  below  the  coincidence  frequency.  The  directivity 
function  is  highly  influenced  by  the  additional  line  stiffness  above 
the  coincidence  frequency  only  near  6 = 90° , where  the  coincidence 
angle  decreased  by  as  much  as  38°  at  SI  = 1.5  and  by  27° 
at  SI®  2.5.  The  peak  pressure  can  be  reduced  by  as  much  as  6 dB  at 
SI  ■ 1.5  and  11  dB  at  SI  ■ 2.5. 

The  radiated  power  from  a plate  with  an  attached  line 
stiffness  reduces  by  only  2.5  dB  below  coincidence  and  by  4 dB  above 
coincidence . 

(iii)  Influence  of  an  Attached  Beam:  In  general,  the  attached 

beam  influences  the  radiated  pressure  spectrum  at  <j>  ■ 0°  only  through 
its  mass,  as  was  discussed  in  Item  (ii)  above. 

The  directivity  function  is  slightly  influenced  by  the 
attached  beam  below  the  coincidence  frequency.  The  directivity 
function  is  highly  influenced  by  the  attached  beam  above  the  coincidence 
frequency,  where  the  coincidence  angle  could  decrease  by  as  much  as 
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30°  at  ft  = 1.5  and  21°  at  51  * 2.5.  It  should  be  noted  that  the 
reason  for  the  shift  is  primarily  due  to  the  stiffness  of  the  beam 
(decreasing  the  angle)  than  to  the  mass  of  the  beam  (increasing  the 
angle) . The  peak  pressure  can  be  reduced  by  as  much  as  6 dB  at 
0 * 1.5  and  12  dB  at  «*  2.5. 

The  radiated  power  from  a plate  with  an  attached  beam 
generally  decreases  with  increasing  mass  and  stiffness  of  the  beam. 

The  reduction  in  radiated  power  reaches  5 dB  below  coincidence,  which 
is  primarily  influenced  by  the  mass  of  the  beam  rather  than  its 
stiffness.  The  radiated  power  reduction  reaches  8 dB  above 
coincidence  which  is  primarily  influenced  by  the  stiffness  of  the 
beam  rather  than  its  mass, 

4.2.4  Conclusion.  The  problem  of  reduction  of  the  radiated 
noise  from  plates  by  attaching  vibrating  machinery  to  the  beam  was 
shown  to  be  successful.  This  study  was  concerned  with  investigating 
the  influence  of  the  mass  and/or  stiffness  of  a beam  on  the  radiated 
sound.  Significant  conclusions  can  be  drawn  from  this  study, 
especially  in  regard  to  the  noise  reduction  problem.  The  main 
conclusions  are: 

(i)  For  excitation  frequencies  below  the  coincidence 
frequency,  the  peak  pressure  occurs  at  $ = 03,  which  can  be  reduced 
significantly  (approximately  8 dB)  by  using  a massive  beam  with  just 
the  necessary  stiffness  dictated  by  structural  design  considerations. 
The  radiated  power  was  similarly  reduced  by  this  mechanism- 

(ii)  For  excitation  frequencies  above  the  coincidence 
frequency,  the  peak  pressure  occurs  at  a new  coincidence  angle  <f^. 

The  peak  pressure  can  be  reduced  significantly  (approximately  12  dB) 
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by  Increasing  the  stiffness  of  the  beam.  The  angle  shift  is  not 
important,  considering  the  large  decrease  in  the  peak  pressure 
accompanying  this  shift.  The  reduction  in  radiated  power  can  be 
reduced  by  as  much  as  8 dB  by  the  use  of  a moderately  stiff  beam. 

4 . 3 Par-Field  Acoustic  Radiation  from  Acoustically-Faired 

Beam-Reinforced  Plate 

The  beam-reinforced  plate  discussed  in  Section  4.2  is  now 

being  acoustically  faired  to  further  reduce  the  radiated  sound. 

The  particular  method  of  fairing  employed  in  this  study  is  done  by 

attaching  an  infinite  line  sprung  mass  to  each  side  of  the  reinforcing 

beam.  Each  sprung  mass  consists  of  a uniformly  distributed  mass  per 

linear  length  m^  and  spring  constant  per  unit  length  and  located 

at  a distance  d from  the  beam.  Thus,  a parametric  study  of  the 

radiated  sound  from  an  arjustically  faired  beam-reinfdrced  plate  can 

. 2 

be  made  with  the  parameters  defined  by  m^/p^h  representing  the 

2 

non-dimensional  sprung  mass,  w, /w  (03  = k, /m. ) the  nondimensional 

1 o i 11 

natural  frequency  of  the  spring-mass  system  and  k^d  the  nondimensional 
distance  to  the  beam. 

The  expression  for  the  radiated  sound  from  an  acoustically 

faired  beam- reinforced  plate  is  given  in  Equation  (2.74).  The 

expression  contains  a factor  C which  takes  into  account  the  beam 

ms 

as  well  as  the  spring-ma9s  systems.  Since  the  sprung  masses  are 
symmetrically  located  with  respect  to  the  beam,  they  vibrate  in 
phase  with  each  other  and  give  off  a sound  pressure  similar  to  two 

line  monopoles  separated  by  distance  2d.  The  directivity  of  such  a 
source  term  is  exhibited  by  the  term  cos  (k^d  cos  0 sin  $) . 
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If  the  forcing  frequency  oj  is  not  close  to  the  natural 

frequency  and  if  the  separation  distance  is  small  compared  to  a 

wavelength,  the  two  monopole  line  sources  merge  into  one  and  the 

cosine  term  can  be  approximated  by  unity.  The  integrals  I ° and 

I? ° are  approximately  equal  to  I 0 if  k^d  <<  1.  Thus,  the 

effect  of  the  sprung  masses  on  the  directivity  function  is  very 

minimal  when  k^d  <<1.  If  the  distance  between  the  sprung  masses  is 

much  larger  than  a wavelength,  i.e.,  l^d  >>  1,  the  integrals  I ° 

and  I become  negligible  and  the  cosine  function  oscillates 

rapidly.  Thus,  the  effect  of  the  sprung  masses  on  the  directivity 

function  is  again  minimal.  It  is  expected  that,  if  k^d  >>  1,  the 

sprung  masses  are  too  far  to  influence  the  near-field  of  the  force 

and,  hence,  the  radiated  pressure  field  must  be  that  of  a vibrating 

<beam-reinf orced  plate.  This  can  be  shown  analytically  if  one  neglects 

1 2 and  I^°  with  respect  to  I^°  in  Equation  (2.74),  where  the 

expression  for  C becomes  C of  Equation  (3.1). 
ms  ms 

If  the  excitation  frequency  uj  is  close  to  the  tuning  frequency 

u^,  the  expression  for  the  radiated  sound  is  given  by  Equation  (3.5). 

At  such  a frequency,  the  line  sprung  masses  become  simple  supports 

for  the  plate,  allowing  complete  freedom  of  rotation  but  no  transverse 

motion.  Thus,  the  Influence  of  sprung  masses  is  very  strong.  If 

k^d  <<  1,  then  the  sprung  masses,  acting  as  simple  supports,  approach 

the  excitation  point  and,  consequently,  the  radiated  sound  becomes 

negligible.  If  k d >>  1,  the  supports  are  too  far  from  the 

excitation  point  to  influence  the  radiated  field.  The  expression  for 

C in  Equation  (3.5)  reduces  to  C for  an  excited  beam-reinforced 
ms  ms 
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plate  given  in  Equation  (3.1)  if  one  lets  12°  and  I^°  vanish.  It  is 
expected  that  when  k^d  'v  0(1)  the  radiated  power  will  increase 
because  of  the  scattering  of  elastic  waves  at  the  supports. 

Due  to  complexity  of  the  acoustic  fairing  problem,  and  the 
large  number  of  independent  parameters  that  need  varying,  the  geometry 
of  the  beam- reinforced  plate  was  made  constant  The  beam  parameters 
used  in  the  numerical  computations  are  q^  = 1 and  y = 2.7. 

In  the  subsequent  discussions,  numerical  results  for  the 
acoustically  faired  beam-plate  were  compared  to  those  for  the  beam- 
plate  only. 

The  influence  of  the  three  parameters  can  be  best  discussed 

in  the  two  frequency  ranges:  below  and  above  the  coincidence  frequency. 

4.3.1  Excitation  Frequency  Below  the  Coincidence  Frequency 

2 

(ft  < 1).  The  influence  of  the  parameters  u>./u)  . m.  Id  h and  k d 
1 o 1 o o 

on  the  radiated  sound  will  be  discussed  in  relation  to  oj/co^,  where 

is  the  tuning  frequency  of  the  spring-mass  system.  In  the 

frequency  range  below  ft  = 1.0,  the  sprung  mass  located  at  k^d  *>  n/5 

is  always  within,  those  at  k^d  = tt  and  3tt  are  always  outside  a half- 

structural  wavelength  from  the  point  of  excitation. 

(i)  ft  < _<  1:  In  this  frequency  range,  the  excitation 

frequency  is  below  the  tuning  frequency  , where  the  action  of  the 

sprung  mass  is  mass  controlled  and  the  force  transmitted  to  the  plate 

2 

is  proportional  to  m^co  . 

The  influence  of  the  separation  distance  k d on  the  radiated 

pressure  spectrum  at  <p  ■ 0°  is  exhibited  in  Figures  3.13  and  3.14. 

It  can  be  seen  that  a modest  reduction  is  attainable  when  k d * n/5, 

o 

but  actually  increases  when  k d * IT  or  3n . The  influence  of  the 

o 
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separation  distance  k^d  on  the  directivity  function,  shown  in  Figure 

3.2*4,  is  predictably  very  marginal.  The  Influence  of  the  separation 

distance  k^d  on  the  radiated  power  is  shown  in  Figures  3.41  and 

3.42.  The  radiated  power  is  somewhat  reduced  (1  to  2 dB)  when  k^d  = tt / 5, 

but  raised  even  bv  a smaller  amount  when  k d = tt  or  3tt. 

o 

2 

The  influence  of  the  sprung  mass  m^/p^h  on  the  radiated 

pressure  spectrum  at  <p  = 0°  is  shown  in  Figures  3.18  and  3.20. 

2 

Increasing  the  sprung  mass  (m^/p^h  = 1 to  5)  generally  decreases  the 

pressure  when  k^d  = tt/5  (Figure  3.18),  the  reduction  reaches  3 dB  when 

the  sprung  masses  are  close  to  the  beam.  The  sprung  masses  have 

little  influence  on  the  pressure  when  their  distance  is  far  from  the 

beam,  as  is  shown  in  Figures  3.19  and  3.20.  The  influence  of  the  sprung 

mass  on  the  directivity  function  is  minimal,  as  is  shown  in  Figures 

3.27  and  3.28,  The  Influence  of  the  sprung  mass  on  the  radiated 

power  is  shown  in  Figures  3.46  and  3.47.  The  radiated  power  decreases 

with  increasing  sprung  mass  (Figure  3.46''  by  as  much  as  4 dB  when 

the  sprung  mass  is  close  to  the  beam  (k  d = tt / 5 ) . No  noticable 

change  in  the  radiated  power  is  observed  when  the  sprung  masses  are 

increased  and  their  distance  from  the  beam  is  far  (k  d = 3tt)  as  shown 

o 

in  Figure  3.47. 

The  influence  of  the  stiffness  k^  is  indicated  by  the  value 

of  the  tuning  frequency  The  influence  of  the  stiffness  on  the 

pressure  spectrum  at  $ • 0°  is  shown  in  Figures  3.15  and  3.17.  The 

pressure  decreases  by  3 dB  with  increasing  stiffness  for  k d * tt/5, 

while  it  increases  by  3 dB  for  k d ■ n or  3tt.  The  influence  of  the 

o 

spring  stiffness  on  the  directivity  function  is  minimal  as  shown  in 
Figures  3.25  and  3.26.  The  influence  of  the  spring  stiffness  on  the 
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radiated  power  is  shown  in  Figures  3. A3  to  3.A5.  In  general, 

increasing  stiffness  decreases  the  radiated  power  by  2 dB  for  k d * tt/5 , 

while  no  noticable  reduction  was  observed  when  k d = TTor  3tt  , as 

o 

anticipated  earlier  in  Section  A. 3 

(ii)  ft  ^ io^/u 1:  When  the  excitation  frequency  oj  is 

close  to  the  natural  frequency  tu  ^ , the  expression  for  radiated 
sound  is  given  by  Equation  (3.5).  The  discussion  in  Section  A. 3 
concluded  that  if  kd  « 1,  the  radiated  sound  is  negligible  and 
if  k d » 1,  the  radiated  sound  is  independent  of  the  sprung  mass 
system. 

The  influence  of  the  separation  distance  k d on  the  radiated 

o 

sound  spectrum  at  <J>  = 0°  when  w ^ is  shown  in  Figure  3.13.  It  can 
be  seen  that  the  reduction  of  pressure  for  k d = tt/5  reaches  15 
dB  at  ft  * w/cJj  = 1.0.  The  pressure  reduction  is  broadband,  where  at 
least  a 6 dB  reduction  is  achieved  for  frequencies  0.5  < ft  < 1.2. 

The  radiated  pressure  actually  increases  when  k d * Tior  3tt,  where  a 

6 dB  increase  is  noticable.  The  influence  of  k d on  the  directivity 

o 

function  is  shown  in  Figure  3.21.  While  the  directivity  function  is 
not  affected  at  k^d  = tt/5,  except  for  a reduction  in  amplitude,  the 
directivity  functions  for  k^d  = tt  or  3tt  have  new  peaks.  These  peaks 
in  the  azimuthal  plane  of  the  directivity  function  can  be  traced  to 
the  minor  lobes  of  the  directivity  function  of  the  cosine  term. 

The  influence  of  the  separation  distance  k^d  on  the  radiated  power 
near  w 't  is  exhibited  in  Figure  3.A1.  The  radiated  power  was 
reduced  by  as  much  as  20  dB  when  k^d  - tt/5,  and  ft  • w/w^  ■>  1.0. 

The  reduction  is  broadband,  where  3 dB  reduction  in  power  can  be 
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achieved  for  frequencies,  0.5  < fi  < 1.2.  The  radiated  power  actually 

increases  (1  jB)  when  k d = it  or  3tt. 

o 

Since  the  spring-mass  systems  act  as  rigid  supports  when 

w = , they  do  not  influence  the  radiated  sound  at  all. 

(iii)  < ft  5.  1*  In  this  frequency  range,  the  excitation 

frequency  w is  above  the  tuning  frequency  , where  the  action  of  the 

sprung  mass  is  stiffness  controlled,  and  the  force  transmitted  to  the 

plate  is  proportional  to  displacement  at  x = + d. 

The  influence  of  the  separation  distance  k^d  on  the  radiated 

pressure  spectrum  at  <P  = 0°  is  shown  in  Figure  3.12.  In  this  case, 

the  radiated  sound  is  reduced  appreciably  by  as  much  as  10  dB  when 

k^d  = 71,  but  the  radiated  sound  spectrum  increased  by  as  much  as  10  dB 

when  k d = tt/5  or  3tt.  The  influence  of  k d on  the  radiated  power  is 
o o 

minimal  in  this  frequency  range,  as  is  shown  in  Figure  3.40. 

The  influence  of  the  sprung  mass  on  the  radiated  pressure  in 

this  frequency  range  should  be  minimal,  since  the  action  of  the  sprung 

mass  is  stiffness  controlled.  This  was  borne  out  by  the  numerical 

2 

computations  for  to.  /w  =0.1  and  m,  /p  h = 1 to  5 at  k d = tt/5,  tt  and 
1 o 1 o o 

3tt.  These  results  are  not  presented  in  this  thesis  for  lack  of 
interesting  phenomena  and  brevity, and  similar  calculations  on  the 
radiated  power  were  computed  and  found  not  to  vary  when  m^  was 
increased  fivefold. 

The  influence  of  the  stiffness  k^  on  the  radiated  sound 
spectrum  at  $ * 0s  is  exhibited  in  Figures  3.15  to  3.17.  For 
to  increase  from  0.1  to  1.0,  the  stiffness  k^  increased  by  a factor  100. 
The  pressure  spectrum  increased  by  as  much  as  10  dB  when  k d = 7T,  but 
decreased  by  as  much  as  8 dB  when  k d = tt/5  or  3tt.  The  influence 
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of  the  stiffness  k^  on  the  radiated  power  is  shown  in  Figures  3.43  to 
3.45,  for  * 0.1.  In  the  frequency  range  0.1  < £2  < 1,  the 

stiffness  seems  to  have  little  influence  on  the  radiated  power. 

4.3.2  Excitation  Frequency  Above  the  Coincidence  Frequency 
(£2  _>  1):  The  influence  of  the  three  parameters  of  the  sprung-mass 

on  the  radiated  sound  will  be  discussed  in  relation  to  u/U)^  again 

(i)  (a^/w  > £2  >_  1:  In  this  frequency  range,  the  sprung 

mass  action  is  mass  controlled  and  the  transmitted  force  is 
proportional  to  the  inertia  of  the  sprung  mass. 

The  influence  of  the  separation  distance  kQd  on  the  radiated 
pressure  spectrum  at  $ * 0°  is  shown  in  Figure  3.14.  For  £2  < 2.5,  the 
separation  distance  kQd  = tt/5  is  within  a half  structural  wavelength 
from  the  point  of  excitation  in  the  plate.  Thus,  the  radiated 
pressure  spectrum  is  reduced  by  as  much  as  10  dB  for  k^d  « tt/5.  The 
influence  of  the  separation  distance  on  the  directivity  function  is 
exhibited  in  Figure  3.30.  The  coincidence  angle  did  not  change  in  the 
plane  6 * 0°  as  expected,  and  there  is  a reduction  in  the  peak 
pressure  by  4 dB  when  k^d  = tt/5,  but  no  major  change  when  k^d  = tt  or 
3tt.  A major  change  in  the  coincidence  angle  and  the  peak  pressure 
occurred  at  the  plane  0 = 90°.  The  coincidence  angle  shifted  up  by 
10°  for  kQd  = tt/5  and  by  14°  for  k^d  = tt,  both  of  which  are  within 
half  a structural  wavelepgth  from  the  beam.  The  cause  of  the  increase 
is  the  same  as  that  of  an  attached  mass,  as  discussed  in  Section  4.2.2. 
When  the  mass  is  outside  a wavelength,  as  is  the  case  for  k^d  = 3tt,  one 
expects  no  change  in  the  coincidence  angle,  as  is  shown  in  Figure  3.30. 
The  influence  of  k^d  on  the  radiated  power  is  exhibited  in  Figure 
3.42.  The  radiated  power  decreased  when  k^d  = tt/5  by  as  much  as  4 dB 
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for  ft  <_  2.  The  radiated  power  for  k^d  = tt  and  3tt  at  ft  = 1.5  does  not 
change  noticeably. 

The  influence  of  the  mass  in  the  radiated  pressure  spectrum 

at  4>  = 0°  is  similar  to  the  case  when  ft  < 1.  Similar  reductions  were 

attained  (approximately  4 dB)  when  the  sprung  mass  increased  five 

times  at  ft  - 15,  tv, /t*J  =2.5  and  k d = tt/5.  The  influence  of  the 
1 o o 

sprung  mass  on  the  directivity  function  is  minimal  as  to  the  shift 

in  the  angle  of  coincidence  as  well  as  the  peak  pressure.  These 

numerical  results  were  computed,  but  are  not  presented  in  the  thesis. 

Increasing  the  sprung  mass  reduces  radiated  power  at  ft  = 1.5, 

tv^/tv^  = 2.5  and  k^d  = tt/5  by  as  much  as  4 dB,  while  no  such  reduction 

is  attained  when  k d = tt  or  3tt.  These  numerical  results  are  not 
o 

presented  in  the  thesis. 

The  influence  of  the  stiffness  k^  on  the  radiated  pressure 

spectrum  at  <f>  = 0°  is  shown  in  Figures  3.15  to  3.17.  The  radiated 

pressure  reduced  by  14  dB  at  ft  = 1.5  and  2.0  when  k d = TT/5  and 

o 

the  stiffness  increased  six  fold.  The  pressure  decreased  by  only  2 

dB  when  k d = tt  and  3n  for  the  same  increase  in  the  stiffness.  The 

o 

influence  of  the  stiffness  is  exhibited  in  Figures  3.31  and  3.33. 

In  both  figures,  one  can  see  that  the  coincidence  angle  did  not 
shift  at  0 =0°.  The  coincidence  angle  did  increase  at  0 = 90°, 
due  to  the  additional  mass  carried  by  the  plate,  especially  when  the 
sprung  mass  fell  within  a structural  wavelength  (k^d  = tt/5. 

Figure  3.31)  while  it  did  not  change  much  when  k^d  = 3tt,  see  Figure 
3.33.  The  influence  of  the  stiffness  k^  on  the  radiated  power  is 
shown  in  Figures  3.43  to  3.45.  The  radiated  power  decreased  by  as 
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much  as  4 dB  at  ft  = 1.5  when  k^d  = tt/5  but  changed  very  little  when 

k d » 7T  or  3n. 
o 

(il)  ^ ft  1:  When  the  excitation  frequency  approaches 

the  tuning  frequency  each  sprung  mass  becomes  a simple  line 

support  for  the  plate.  The  discussion  follows  the  same  lines  as  the 

one  given  in  Section  4.3.1.  The  radiated  pressure  is  then 

independent  of  the  mass  or  stiffness  of  the  system  but  is  dependent 

on  the  value  of  k d.  If  k d <<  1,  the  radiated  sound  becomes 
o o 

negligible,  and  if  k^d  >>  1,  even  the  distance  kQd  does  not  influence 
the  radiated  sound. 

The  influence  of  k d on  the  radiated  pressure  spectrum  at 

o r 

4)  = 0°  when  u)  ^ U)^  is  exhibited  in  Figures  3.13  and  3.14.  In  each 
case,  the  radiated  pressure  was  reduced  considerably  when  k^d  = tt/5, 
by  as  much  as  15  and  22  dB  at  * ft  = 1.0  and  2.5,  respectively. 

The  radiated  pressure  increased  when  k^d  = 3tt,  where  the  system  is 
located  well  outside  a structural  wavelength  for  all  frequencies 
up  to  tl  s 2.5,  by  as  much  as  8 dB . 

The  influence  of  the  separation  distance  k d on  the  directivity 
function  is  exhibited  in  Figures  3.22  and  3.23.  The  coincidence 
angle  and  peak  in  the  pressure  is  completely  absent  from  the  case 
kQd  * tt/5  at  6 * 0°  and  9 = 90°  and  the  frequencies  ft  = 1.5  and 
2.5.  The  influence  of  larger  values  of  k^d  is  negligible  when  9 = 0° 
but  can  shift  up  or  down  depending  on  the  value  of  k^d.  One  should 
note  that,  unlike  the  discussion  in  Item  (i)  above,  the  sprung 
mass  no  longer  contributes  to  the  surface  mass  density  of  the  plate 
and,  hence,  a physical  reasoning  cannot  be  made  to  account  for  the 
shift  in  the  angle  of  coincidence.  The  peak  pressure  reduced  by  as 
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much  as  18  dB  for  0=0°  and  25  dB  for  G = 90°  when  k d = t/5 

o 

and  ft  = 1.5  or  2.5,  while  it  increased  by  1 dB  at  ft  =1.5  and  7 

dB  at  ft  = 1.5  for  k d = 3tt.  The  influence  of  k d on  the  radiated 

o o 

power  is  shown  in  Figures  3.41  and  3.42.  The  radiated  power  generally 

decreased  by  16  dB  at  ft  = 1 and  23  dB  at  ft  = 2.5  when  k d = tt/5 . 

o 

The  radiated  power  actually  increased  when  k d - u or  3ir  by  as  much 
as  3 dB  at  ft  = 1.5  and  2.5. 

(iii)  ft  > ^ ^is  frequency  range,  the  sprung 

masses  are  stiffness  controlled.  As  o j/uj^  increases  without  limit, 
it  can  be  seen  from  Equation  (2.74)  that  the  effect  of  the  sprung 
masses  will  be  negligible  and  the  acoustically  faired  reinforced 
plate  will  behave  as  if  the  sprung  masses  are  absent. 

The  influence  of  the  separation  distance  k^d  on  the  radiated 
spectrum  at  = 0°  is  illustrated  in  Figure  3.13.  Comparison  of  the 
spectra  when  cu^/oj^  = 1-0  shows  that  regardless  of  the  distance  k^d, 
they  all  approach  the  spectrum  for  beam-reinforced  plates  from  above 
as  ft  increases.  The  influence  of  kQd  on  the  directivity  function  is 
shown  in  Figure  3.35.  In  this  case,  when  ft  = 2.5  and  = 1-0,  the 

distance  k^d  has  negligible  effect  on  the  angle  of  coincidence  as 
well  as  the  peak  pressure.  The  influence  of  the  distance  k^d  on 
the  radiated  power  is  shown  in  Figure  3.41.  Again,  the  radiated 
power  is  not  altered  significantly  for  all  values  of  k^d. 

The  influence  of  the  sprung  mass  m^  on  the  radiated  pressure 
spectrum  at  <p  = 0°  becomes  minimal  as  ft  >>  as  can  be  seen  in 

Figures  3.18  through  3.20,  where  the  spectra  for  approach  the 
beam-plate  spectrum  from  above.  The  influence  of  the  mass  m^  on 
the  directivity  function  is  exhibited  in  Figures  3.37  and  3.39. 
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In  each  case,  the  mass  has  no  influence  on  the  angle  of  coincidence 
or  on  the  peak  pressure.  The  influence  of  the  sprung  mass  on  the 
radiated  power  also  diminishes  as  £2  increases,  as  is  shown  in  Figure 
3.46,  where  the  power  spectrum  for  the  sprung  mass  approaches  the 
power  spectrum  for  the  beam-plate  from  above. 

The  influence  of  the  stiffness  on  the  radiated  pressure 
spectrum  at  $ = 0°  also  becomes  negligible  as  £2  becomes  larger  then 
as  is  shown  in  Figures  3.16  and  3.17.  The  influence  of  the 
stiffness  on  the  directivity  function  is  exhibited  by  Figures 
3.36  and  3.38.  The  directivity  function  for  = 1.0  is  very 

close  to  that  for  = 0 (no  sprung  mass),  where  the  angle  of 

coincidence  and  the  peak  pressure  changed  very  little.  The  influence 
of  the  stiffness  on  the  radiated  power  is  shown  in  Figures  3.44  and 
3.4 j where  the  comparison  is  to  be  made  on  the  power  spectrum  for 

* 1.9.  In  each  case,  the  power  spectrum  approaches  that  of  a 
beam- reinforced  plate  from  above  as  £2  increases. 

4.3.3  Summary . In  this  section,  a summary  of  the  discussions 
carried  out  in  Section  4.3  is  presented.  The  summary  will  be  present 
in  the  three  frequency  ranges  <1,  ^ 1,  and  > 1. 

(i)  > £2:  In  this  frequency  range,  the  sprung  mass's 

action  is  mass  like. 

Below  the  coincidence  frequency,  the  most  important 

parameter  that  influences  the  radiated  sound  is  the  separation 

distance  k d.  If  k d • TT/5  is  well  within  half  a structural 
o o 

wavelength  for  £2  > 0.1,  the  radiated  pressure  and  power  spectrum 
decreases  with  increasing  mass  m^  or  stiffness  k^.  If  k^d  « 3tt  , 
well  outside  a structural  wavelength  for  £2  < 2.5,  the  radiated  sound 
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usually  increases  with  increasing  sprung  mass  m^  or  stiffness  k^. 

The  influence  of  all  the  parameters  on  the  directivity  function  is 
minimal . 

Above  the  coincidence  frequency,  the  same  can  be  said  for  the 
radiated  sound  pressure  and  power  spectrum.  The  only  change  that  is 
different  is  the  influence  of  the  sprung  mass  on  the  angle  of 
coincidence.  The  angle  usually  increases  if  kQd  is  within  a 
structural  wavelength  from  the  beam,  primarily  due  to  the  added  mass 
effect,  and  is  negligible  if  k^d  is  outside  a structural  wavelength. 

(ii)  ^ &'•  If  the  excitation  frequency  w is  close 

to  the  natural  frequency  of  the  sprung  mass  co^,  the  sprung  masses 
become  simple  supports.  Hence,  the  mass  m^  and  stiffness  k^  has 
no  influence  on  the  radiated  pressure.  Thus,  if  the  distance  k^d  is 
small  compared  to  structural  wavelength,  then  the  radiated  pressure 
spectrum  can  decrease  by  as  much  as  15  dB  below  coincidence  and  by 
22  dB  above  coincidence.  The  reduction  is  broadband  about  the 
frequency  u>  * , extending  to  over  an  octave  for  a 6 dB  reduction. 

The  radiated  pressure  actually  increases  by  as  much  as  6 dB  if  k^d 
represents  a distance  outside  a half  structural  wavelength.  The 
reduction  in  the  pressure  is  achieved  over  the  whole  azimuthal 
plane,  so  that  the  radiated  power  does  reduce  by  as  much  as  20  dB  below 
coincidence  and  16  dB  above  coincidence. 

(iii)  < this  frequency  range,  the  sprung  mass’s 

action  is  stiffness  controlled.  As  the  frequency  of  excitation 
increases  over  the  tuning  frequency  the  influence  of  the 

sprung  mass  diminishes,  as  was  expected  when  the  analytical  solution 
was  examined.  Thus,  the  radiated  pressure  and  power  spectra  for  an 
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acoustically  faired  beam-reinforced  plate  approach  these  for  a 
beam-reinforced  plate  from  above  as  ft  >>  . The  directivity 

function  is  also  unaffected  when  ft  >>  cu,  /to  . 

1 o 

4.3.4  Conclusion . The  problem  of  reduction  of  the  radiated 
noise  from  reinforced  structures  was  attempted  by  employing  a method 
of  acoustic  fairing  where  an  infinite  line  sprung  mass  is  added  to  each 
side  of  the  reinforcement.  After  an  exhaustive  number  of  parametric 
numerical  results  were  computed,  only  a small  portion  of  which  were 
presented  in  this  thesis,  the  main  conclusions  are  presented  below: 

(i)  To  achieve  a reduction  in  the  radiated  pressure  and 
power  spectra,  it  is  recommended  that  the  acoustic  fairing  is  applied 
to  within  half  a structural  wavelength  from  the  reinforcing  beam. 

(ii)  To  achieve  a broadband  noise  reduction,  the  sprung-mass 
system  employed  in  the  acoustic  fairing  must  be  tuned,  such  that  its 
natural  frequency  is  close  to  the  frequency  around  which  a reduction 
is  desired.  This  tuning  provides  as  much  as  25  dB  noise  reduction 
at  the  desired  frequency.  The  tuning  normally  provides  for  almost 
an  octave,  (centered)  at  the  desired  frequency,  where  at  least 
6 dB  noise  reduction  is  achieved. 

(iii)  Structural  discontinuities  of  any  kind  usually  increase 
the  radiated  pressure  and  power  in  the  frequency  range  below  ft  = 2.5 
if  they  are  located  outside  a structural  half  wavelength.  Thus, 
such  discontinuities  are  to  be  avoided  if  at  all  possible,  if  one 
desires  to  reduce  the  radiated  noise  below  ft  « 2.5. 
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ii  Supplementary  notes 
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12  SPONSORING  MIL!  T AR  y ACTIVITY 

Naval  Ordnance  Systems  Command 
Department  of  the  Navy 

*The  radiated  sound  from  a beam  reinforced  plate  excited  by  a point  force  is 
investigated.  The  point  force  is  located  on  the  reinforcing  beam,  simulating  a standar 
method  of  supporting  vibrating  machinery  by  plate  structures.  The  radiated  pressure 
spectrum  of  excited  reinforced  plates  was  reduced  by  as  much  as  6 dB  for  moderately  stifj 
beams.  The  radiated  power  generally  decreased  with  increasing  mass  and  stiffness  of  th<j 
reinforced  beam,  reaching  by  as  much  as  8 dB  for  excitation  frequencies  above  the  coin- 
cidence frequency.  The  angle  of  coincidence  usually  increased  with  increasing  mass  and 
decreased  with  increasing  stiffness  of  the  reinforcing  beams.  The  peak  pressure  at  the 
coincidence  angle  could  be  reduced  by  as  much  as  12  dB  by  judicious  choice  of  the  geom- 
etry of  a reinforcing  beam,  whose  dimensions  are  generally  dictated  by  structural 
design  considerations. 

A method  of  adding  other  structural  discontinuities  to  further  reduce  the  radiated 
noise  was  also  investigated.  This  is  referred  to  as  Acoustic  Fairing  in  this  thesis. 

The  procedure  thus  employed  consists  of  a pair  of  infinite  line  sprung  masses,  where 
one  if  located  on  each  side  of  the  reinforcing  beam.  The  radiated  sound  from  point 
cxc  in  I reinforced  plates  with  acoustic  fairing  was  investigated  and  a parametric  study 
was  mpleted.  It  was  found  that  a modest  2-4  dB  reduction  in  the  radiated  pressure 
and  p wer  can  T>e  attained  if  the  sprung  mass  is  located  within  a half  structural  wave- 
i«*ngt!  from  the  beam.  In  addition,  if  the  tuning  frequency  of  the  sprung  mass  is  close 
to  f-  excitation  frequency,  reduction  of  the  radiated  pressure  and  power  of  the  order 
IS-.'.'  dB  were  realized.  The  reduction  thus  achieved  was  also  broad  band,  extending  oveij 
an  octave  centered  at  the  tuning  frequency.  If  the  sprung  mass  is  located  outside  a 
structural  wavelength  from  the  beam,  an  increase  in  the  radiated  pressure  and  power 
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13.  ABSTRACT  (continued): 

results,  primarily  through  scattering  of  elastic  waves  at  the  sprung 
mass  junctions.  The  effectiveness  of  the  sprung  masses,  regardless 
of  their  location  with  respect  to  the  beam,  is  vastly  reduced  when 
the  excitation  frequency  is  much  higher  than  the  tuning  frequency. 
The  effectiveness  of  the  sprung  masses  also  diminishes  when  they 
are  located  too  many  wavelengths  away  from  the  beam. 
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